5. ESTRUTURAÇÃO DO MODELO 


Cada unidade de exploração terá de ser desmontada de baixo para cima em sucessivos «cortes» ou 
«plantas», de altura fixa ao longo de toda a vida do jazigo e comum em todas as suas unidades de explo- 
ração, fixação essa que é consequência da normalização do material de perfuração a usar e da forma de 
carregamento dos tiros de cada pega de fogo. A tonelagem assim abatida, em cada «corte» ou «planta» do 
«cut and fills, é função apenas da secção do «céu» de cada câmara de desmonte e não da secção do «céu» 
da unidade de exploração e será constante em cada«corte» ou «planta» dessa mesma câmara, pois em ter 
mos médios genéricos pode-se dizer que a secção de cada câmara é constante ao longo de toda a sua vida 
(a excepção a esta generalidade são as «franjas» do jazigo). 

Por motivos de segurança é conveniente ter-se sempre cheias de minério as toldas de descarga, isto é, 
ao dar-se início ao carregamento do minério abatido em dado «corte» ou «plantas, a tolda da câmara 
deverá estar ainda cheia com o último minério carregado do «corte» ou «planta» anterior e, ao terminar 
o carregamento do minério abatido, deve ficar por retirar da tolda o minério suficiente para que a tolda fique 
cheia até ao início do carregamento do «corte» ou «planta» seguinte. Assim, desprezando o início e o final 
da vida de cada piso, é lícito considerar que em cada fase de carregamento de minério, a quantidade total 
«carregada» da unidade de exploração é exactamente a quantidade abatida por «cortes ou «planta», para 
além apenas do acréscimo de altura das toldas de «corte» para «cortes e que é verdadeiramente insigni- 
ficante em face da tonelagem movimentada em cada fase de carregamento. 

Portanto, a cada fase de carregamento na câmara corresponde igual período de carregamento de 
vagonas na galeria de rolagem, vagonas essas que serão comboiadas para o poço de extracção. A este 

T, 
poço chegará, com origem na unidade X,, a tonelagem ue em cada fase de carregamento de cada uma das 
câmaras da unidade, admitindo igualdade no valor da secção de cada uma daquelas duas câmaras de des- 
monte (a excepção a esta generalidade são ainda as «franjas» do jazigo). 

Adicionalmente é condição desejável a utilização de frota de rolagem constante, quer de locomotivas, 
quer de vagonetas, em cada piso, mas compatível com a condição imposta da constância de e, com plena 
utilização durante toda a vida do piso (ou parte do mesmo). 

Nestas condições, é aplicável o teorema: 


«As condições de exaustão do jazigo (piso, no nosso caso) e a de manutenção da frota constante 
implicam que o centro de massas do jazigo se mantenha constante», 


Este teorema estã devidamente demonstrado pelos autores que temos vindo a citar na sua 2.º comu- 
nicação ao I CHILAGE, intitulada «Um método de programação de desmontes>», e é aplicável ao caso que 
apontamos desde que se admita o que se tem vindo a expôr. 

Mas a estacionaridade da frota de rolagem utilizada não deve ser o único objectivo, pois interessa 
que o desmonte seja conduzido por forma a manter estacionários outros parâmetros como, por exemplo, o 
teor de um dado elemento no tal-qual extraído do piso. E os mesmos autores, no mesmo trabalho citado, gene- 
ralizaram a aplicação do teorema anterior para quaisquer outros parâmetros, através da consideração de 
que o valor médio do parâmetro escolhido (qualquer que seja o parâmetro, repita-se), com o minério todo 
do piso explorado por um único acto de exploração, define um a como a distância a uma origem arbitrária 
de um novo centro de massas no espaço dos valores do parâmetro escolhido (teores de um dado elemento, 
por exemplo); portanto, para que o centro de massas, nesse espaço se mantenha estacionário (isto é, após 
cada acto de exploração, o piso quanto ao parâmetro escolhido, seja uma réplica do piso original) e para 
que o valor do parâmetro escolhido do minério trans portado se mantenha constantemente igual a a, é neces- 
sário que, em cada acto de exploração se verifique: 


e] 


ze, 


E assim se construirá um sistema de equações que representa o conjunto de relações estruturantes 
que definem as associações permitidas entre unidades de exploração de cada piso, reduzindo substancial- 
mente o número de combinatórias possíveis. Está assim ultrapassado o método aproximativo clássico, que 
supõe a massa do piso concentrada no seu centro de massas e arbitra a distância desse ponto ao ponto de 
descarga (poço de extracção), resolvendo, em primeira aproximação, o problema da definição da frota, mas 
nada dizendo sobre a distribuição da frota de modo a manter constante o centro de massas do piso nem 
as características médias do material desmontado e extraído; para além de passarmos a possuir um método 
científico, abandonando o carácter aproximativo do método clássico. Note-se ainda que, quanto maior for 
o número de parâmetros considerados, maior é a complexidade do cálculo dos fluxos de minério das diver- 
sas unidades do piso, e a aproximação clássica é incapaz de resolver o problema de forma aceitável. 
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Por outro lado, o método científico exposto permite sempre, pelo menos teoricamente, encontrar a 
solução ideal para aquele problema, bastando para tal resolver e encontrar a solução do sistema de equa- 
ções constituído pelas relações estruturantes que definem as associações permitidas entre unidades de explo- 
ração de cada piso, nas condições já referidas anteriormente, registe-se que este será um sistema de q 
equações a p incógnitas, sendo 5 o número de parâmetros considerados. Teoricamente haverá solução para 
este sistema de equações, mas é evidente que a sua compatibilidade prática, em geral, não existirá, pelo 
que a resolução de um problema concreto ter-se-á de fazer pela seguinte forma: 


a) à custa de uma das equações apenas, definir as combinações possíveis, 


b) das combinações possíveis anteriores considerar apenas aquelas que satisfaçam a uma das equa- 
ções restantes, obrigatoriamente mais fraca que a equação considerada em a); 


c) e assim sucessivamente até ao ponto desejado, mas tendo em atenção que, quanto mais avan- 
carmos na restrição da combinatória, menor será o número de combinações possíveis. 

Assim, a sensibilidade do calculador manifestar-se-ã na escolha do ponto de paragem na sucessiva 
restrição da combinatória, mas haverá sempre a possibilidade de se recorrer ao aumento do desvio admissi- 
vel (1) para cada parâmetro, de forma a elevar o número de combinações utilizáveis. 

Para além das restrições assim impostas, há sempre ainda a possibilidade de impôr outras restrições 
de natureza tecnológica, que também irão diminuir o número de combinações permitidas. E o caso, por 
exemplo, da existência obrigatória de uma pá auto-carregadora em cada unidade de exploração e que não 
poderá trabalhar simultaneamente nas duas câmaras da sua unidade de exploração, pelo que serão elimi- 
nadas combinações de unidades que conduzam a excessivas ou deficitárias necessidades de carregamento. 


6. PLANIFICAÇÃO DAS FROTAS E SUA DISTRIBUIÇÃO 


Já anteriormente referimos, mais de uma vez, o interesse e a necessidade de determinação da frota 
ideal para a exploração racional de um jazigo, logo de um piso no nosso caso. E temos de distinguir os 
dois tipos de frota existentes numa lavra subterrânea: a de carregamento e a de rolagem. Ambas impor- 
tantes e com o mesmo grau de indispensabilidade, Vamos referi-las separadamente, embora para ambas o 
problema se possa equacionar com um objectivo idêntico: estabelecer a frota necessária para a taxa e ren- 
dimento constantes, desmontar totalmente o piso em causa e, posteriormente, distribuir por esse piso, ao 
longo do tempo, essa frota, ou seja, estabelecer as taxas de desmonte em cada unidade de exploração de 
forma a manter constantes as características consideradas na estruturação do nosso modelo tal como foi 
feito em 5. 


6.1. FROTA DE CARREGAMENTO 


Como condição prévia aceitamos que a frota é constituída, em cada piso, qualquer que seja o seu 
número de unidades de exploração, por pás auto-carre gadoras iguais. 

Segundo um esquema do tipo do constante da fig. 5, um ciclo da pá auto-carregadora, na fase de 
carregamento de minério, na unidade X, do piso i, é definido por: 


RV + Fm 


x V o V 
e a capacidade de transporte, na unidade de tempo, dessa pá auto-carregadora é dada por: 


EO as gnv 


i E. 
£ ly Ea V = a 


(1) — O desvio y admissível, por exemplo, para o caso do parâmetro a é introduzido na equação definidora de a 
do seguinte modo: 


be 
DO, 


“o 
=e, 
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e o número de pás auto-carregadoras é dado por: 


— na unidade X,: 
“ SM tg * 6 ee bh 
á E 4 x x 
Dj e a = = 
1 Pi 7 “Av 
x 
-——no piso i: 
= V e Tê 
4 nim = 2% 
x E 1] E xá 


mas este somatório ba & “, é reduzido, no máximo, a dois termos, porque, para se manter a produção 
com características médias semelhantes ao longo do tempo deve bastar trabalharmos com o mínimo de 


P- ponto de carga; O— tolda (ponto de des- 
carga ) 
OP- distância a percorrer em vazio 
PO- distância a percorrer em carga 
| = OP + PO —distância de carregamento a 
(CAMARA A) considerar na unidade X, 


Fig. 5 


Esquema de um ciclo de uma pá auto-carregadora 
na fase de carregamento na câmara A da uni- 
dade X, 


duas unidades de exploração diferentes simultaneamente (o uso de maior número de unidades deverá compli- 
car a rolagem no piso, exigindo grande disciplina de utilização do equipamento, sinalização especial, etc...). 


Eliminando Es teremos: 


1 e tvrrz eh. 
Dn; = nm = E 6 tn del) ———— 
nv pmnv 
pois por definição, e, = IJ e, 
A frota constante no piso obriga-nos a que n, = constante, para o que terá que ser E e, 1, = cons- 


tante — Z, constante esta (Z) que irá, verdadeiramente, determinar a frota, pelo que há que lhe atribuir um 
valor para esse efeito. Essa atribuição é feita considerando que a exploração do piso i em questão é uma 
sucessão no tempo de actos de exploração e que em cada um desses actos deve ser cumprida a mesma 
constância de 4; 

Se assim for é sempre possível explorar o piso até à exaustão (|). A totalidade dos actos de explo- 
ração, conduzindo à exaustão do piso, exprime-se por: 


nNE=% Ra la 


(1) — É evidente que não se explorarão os pilares, maciços de protecção, etc., não considerados quando se fala em 
exaustão de um piso e não contabilizados nas reservas exploráveis, 
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onde m representa o número de actos de exploração: 


> Ei, d, + Ei, LE e us Te Ei ho =3 KR. ho. 


mas como já temos: 


É 6, = € 
me =: KR, 
e 
pode-se escrever: 
Z 3 R., =€ 2 RB, ho. 
Z=e L ousejalJe,h,=e Li 
donde: 
oe, h 
L, = - 
Ê;, 
e a frota necessária fica definida por: 
e (vt +L 
i i 
= EP 3 (A) 


| mv 


pois já tinhamos provado que; 


e; tv + É Ci E 
n; = mem - 
nv 


Para esta definição do cálculo da frota necessária (n,) para desmontar todo o piso i ser coerente, 
resta demonstrar a afirmação feita anteriormente de que era possível desmontar a totalidade do piso por 
actos elementares de exploração, obedecendo à condição 


== e, 1, = constante 


ou: 


E Ci hs 
L, = 
É; 


mas tal já foi feito pelos autores que temos vindo a citar na sua comunicação «Um método de progra- 
mação de desmontes» ao I CHILAGE, para o que bastou partirem da definição do centro dz massas de um 
jazigo (piso, no nosso caso) no início do seu desmonte: 

Y 


ad Ro h=R, L onde Y é o número total de unidades de exploração em consideração no piso i. 
X= 


Em resumo, a frota escolhida de acordo com a relação (A), que representa a frota necessária para 
desmontar um piso é a estritamente necessária e trabalhará a rendimento constante, desde que se tomem 
disposições para permitir que o valor de y adoptado em (A) seja o efectivamente verificado na prática; 
os actos de exploração podem incidir, em qualquer momento, sobre quaisquer unidades de exploração desde 
que a condição Z e, 1,=e; L, seja cumprida e, como corolário, L, manter-se-á constante ao longo da vida 
do piso. 

Temos, assim, a frota de carregamento definida, constante ao longo do tempo de vida do piso i, capaz 
Ge desmontar o piso à taxa de exploração e. Resta definir a distribuição dessa frota, mas para tal, intrin- 
secamente, basta que essa distribuição satisfaça a relação: 
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já atrás referida. Mas para que o piso forneça sempre minério com características constantes e iguais à 

média global, haverá que satisfazer outros requisitos. Estez, normalmente, dizem respeito a rendimentos 

em peso e a teores de dados elementos químicos, mas sejam eles quais forem são de considerar cada um 

como um dos parâmtros a largamente debatidos atrás quando tratámos da estruturação do nosso modelo. 
Por exemplo, em matéria de um rendimento em peso 4 como parâmetros teriamos: 


e, de seguida, em matéria de um teor de concentrados € como parâmetro teríamos: 
E 


2& e TA t , 


A x 


ly Miy 


e assim sucessivamente até possuirmos o sistema de equações falado quando da estruturaço do nosso modelo, 
sistema esse que uma das equações será a ligação: 


e recorda-se que para satisfazer simultaneamente todas as equações do sistema deverá ter de se admitir uma 
certa variação em todas ou em algumas das características consideradas (é a questão falada de termos 


t+ : Yo 


6.2. FROTA DE ROLAGEM 


Como condição prévia teremos de aceitar novamente que, em cada piso, o equipamento disponível 
(locomotivas e vagonas) são iguais. 

A fig. 6 procura representar um ciclo de rolagem normal no piso i, ciclo esse que é definido, em relação 
à unidade de exploração genérica X., por: 


R (hk+1)d, Rv t+ (bh, 
DR A x x 
é | ç | Ea + tg : + ly = 
e. x Va V, h v, h v, 
OS A A 
| | || 
A | B Me - [ e Il--ponto de chegada das toldas das duas câma- 
| | | | X; em) 
| 4 | | ras de cada unidade de exploração. 
ne ae aam con Estes mê a ecoa 
PISO | p— IH | UM a 3H III a,, distância de rolagem para o minério da uni- 
E di, oi dade X.. 
POÇO 
a TR ACÇÃO 


Fig. 6 


Esquema de um piso vulgar 


e a capacidade de transporte, na unidade de tempo, de uma locomotiva é dada por: 


Por diversas razões, entre elas de segurança, deve-se procurar que cada piso tenha em funcionamento 
na mesma unidade de tempo apenas uma locomotiva, a qual durante o relevo estará em constante vai-vém, 
podendo completar N. ciclos ou composições no piso i: 


N=2N. +N,. + ENG E ds EE Ma 
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Cada composição deverá, quando em carga, transportar a capacidade máxima de transporte da loco- 
motiva logo, para obtermos na arreação do piso o nosso objectivo e,, terão de lã chegar N, locomotivas: 


' s 
G=N p= No P, 


e porque e, = LZ e, teremos: 
e g Kv,t K + 1jd 
: E; e | e . 1 1 + + ) ii 
N, —— T =— == “ 4 4" — 1 n > É, É | =—— ) 1 z É, K == 
P, Pp”, mm Y x x x Av 
x 

| K Vita + (K + 1) di. | ekvs+Hs “4. (K+1) di ekKvt+(KA+ 1) 26, d 

=" ITD Ds > >—>—>—————— 
di en Kov, ow Ky, qn Koyv, 


Segundo o raciocínio que tem vindo a ser aplicado N, terá de ser constante, o que implica a cons- 
tância de 
* e, d, =E 
e, por analogia com o que se passou com a frota de carregamento, teremos: 


mZ=2:R, d 


| 


onde m continua a representar o número de actos de exploração. 


Desenvolvendo: 


e o número de locomotivas necessárias fica definido por: 


E A Kv tH(K+1) e D e [Kv +(K+1) D] 


(B) 


admitindo-se como verdadeira já a possibilidade de desmonte da totalidade do piso por actos elementares 
de exploração e devendo apresentar-se o mesmo tipo de considerações apresentadas quando abordâmos a 
frota de carregamento (my efectivamente verificado na prática e, consequentemente, teremos D, constante 
ao longo da vida do piso). 

Para este caso da rolagem, resta determinar apenas o número de vagonas que cada composição 
comporta, o que é fácil pois tendo sempre vagonas iguais, de capacidade útil, em peso, q”, teremos ape- 


nas que dividir q' (capacidade da locomotiva) por q', Isto significa que cada composição transportará 


vagonas, havendo que determinar, em função de N, e do tempo de carga e descarga de cada vagona, 
o número de vagonas que existirá no piso, pois se admite que a locomotiva não espera, isto é, chega com 
uma composição carregada à arreação do poço e logo deverá partir com outra composição vazia e vice- 
-versa junto ao local de carregamento das vagonas de forma a respeitar as condições introduzidas (t, ser 
constante, termos apenas uma locomotiva no piso, etc.). 


E evidente que na busca da melhor solução para esta questão se encontra como número de com- 
posições mínimo possível 2+6, sendo O o número de unidades de exploração em actividade de carrega- 
mento simultaneamente, E este número é dito mínimo porque para tal é necessário que seja possível carre- 

t 


gar uma composição durante um ciclo de rolagem menos = admitindo a igualdade de tempo gasto no 
z 


TÉCNICA 451/452 281 


carregamento de uma composição e da sua descarga na arreação, isto é, tanto a carga como a descarga 
terão de ter duração inferior a: 


t 


EL d le k + 1 
+ 1, + ty, =p — E E —.. ade n 
E x x 2 V K v * Y E 
| 1 1 
Este modelo é um dos vários possíveis, pois pode admitir-se, por exemplo, a assistência das loco- 
motivas ao carregamento das composições, entre outros modelos. 


Como nota final, no que diz respeito à frota de rolagem, não queremos deixar de referir o cálculo 
do múmero de vagonas que há que encher para a satisfação da totalidade de minério carregado pelas 
pás auto-carregadoras em dado período de tempo, Em cada ciclo da pá auto-carregadora é colocado, na 


tolda de descarga, minério suficiente para o enchimento de —, vagonas, número de vagonas este que 


' 
à 


a 
pode ou não ser suficiente para o preenchimento de uma composição de vagonas. Como cada composição 
tem + Vagonas teremos que: 


se p>qy' em cada ciclo da pá auto-carregadora encher-se-á uma composição; 
se y>y' em cada ciclo da pá auto-carrega dora encher-se-á mais do que uma composição (exac- 


tamente — composições); e 
E 
— tm <' em cada ciclo da pá auto-carregadora encher-se-ã menos do que uma composição (exacta- 


mente há de uma composição). 

% 

A situação ideal seria a de q =q', pois isso significaria que de cada vez que se completasse um 
ciclo da pá auto-carregadora teriamos de ter carregado integralmente uma composição de vagonas, para 
se respeitar o princípio de segurança de manutenção das toldas de descarga de minério sempre cheias. 


Mesmo matematicamente é esta a hipótese mais cómoda, pois a igualdade de q e q' arrastará a possi- 
bilidade de fácil relacionação entre n, e N, tal como estes valores se expressam em (A) e (B), tudo depen- 
dendo das velocidades v e v, e das distâncias médias L e D. 

A 2º hipótese (> q') não é usual, pois ela significa que as composições de vagonas são muito 
pequenas, o que é anti-econóômico, pois exigiria grande número de locomotivas (todas de pequena capa- 
cidade) e aumento da mão-de-obra respectiva. 


E 


A 3.º hipótese (p<p') é a mais usual e nela verifica-se que o número de ciclos da pá auto-carre- 
gadora para encher uma composição é dado por: 


enquanto que a produção horária de uma composição é: 
pet PD on 
d i E e e E = 
» z E 
X 
desde que =, esteja expresso em minutos. Isto significa que o número de composições necessárias, numa 


1 
hora, para a unidade X,, é dada por: 


r P, O ty N. oO o 


desde que a unidade de tempo h, a que se refere a taxa de desmonte e,., seja a hora. Posto isto, é fácil 


E 
ter o número de vagonas a encher, pois já vimos que ele se relaciona com qw' da seguinte forma: 


2” 
N, =— 
dj 
ct E a E, õ; 
e ' , a É 1 ” E 
então, numa hora, há que encher na unidade X, - DN. Vagonas, a que correspondem SE sea 
60. Nvq | p 60. é 


E 


composições, sendo também estes os números de vagonas e de composições provenientes da unidade x, 
que chegam durante uma hora à arreação do poço de extracção do piso i. 
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7. POÇO DE EXTRACÇÃO 
Seguindo agora raciocinio análogo, teremos como ciclo no poço, para o piso i: 


H, H, n H, TV+2H, 
+ —— sm É 4h É E e 


E = [+ 

Pg Va d 

desde que H, seja a profundidade do piso i e que a velocidade de subida (v, = V) seja igual à veloci- 

dade de descida (v, = V), o que é normal, pois é normal a existência de poços com dois canais tendo skip 

ou uma jaula em cada canal a actuarem em simetria; mesmo que assim não aconteça e exista apenas um 
skip ou uma jaula temos uma velocidade média V válida para aplicação na fórmula. 
Em 4 já tinhamos definido, para o poço de extracção os seguintes rendimentos: 


E, 


-—» (rendimento real, em relação ao piso i, durante o tempo h) 


Ph; — —> (rendimento horário real da extracção, durante o tempo h, em relação ao piso i) 


pelo que a sua capacidade de extracção do piso i na unidade de tempo é dada por: 


= *p 'h o Ta "h; V 
p O ——————= ES 


P; TV+2H, 


Esta capacidade de extracção, numa hora, tem de equivaler à produção horária do piso, produção esta 
proveniente de 6 unidades de exploração em actividade de carregamento simultaneamente. De cada uma 


daquelas unidades X, chegará, numa hora, - composições e a totalidade de composições chegadas à 
arreação do piso, numa hora, será di 

O Ci 

=1 P, 


Como cada composição tem a capacidade útil y' teremos na arreação do piso, numa hora a seguinte 
produção: 


x 


Pp; — p' > 
i=1 PP, 


Para a totalidade do jazigo, com n pisos em actividade servidos por este poço, teremos uma possi- 
bilidade de extracção de: 


2) 
|| 
po t4m 
E 


a qual é repartida pelos diversos pisos em função das características pretendidas para o total extraído, o 
que levou a fixar a distribuição horária da extracção por esses pisos. 

Para cada piso já sabemos as possibilidades de extracção do poco numa hora ou em qualquer outra 
unidade de tempo. 

Não nos parece relevante o prolongamento deste estudo para o poço de extracção, pois ele em qual- 
quer mina bem dimensionada não é o «gargalo» da produção; normalmente esse «gargalo» ou é o equi- 
pamento de carregamento e/ou de rolagem ou é nas instalações a juzante do poço de extracção, 


8. EXEMPLO DE APLICAÇÃO 


O método acima descrito não é ainda utilizado em qualquer das minas onde exercemos a nossa acti- 
vidade profissional actualmente. 

Os resultados práticos da sua aplicação só poderão vir a ser encontrados depois de aplicado o 
método e de os compararmos com a actual metodologia clássica que é seguida, No entanto, métodos 
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semelhantes para a lavra a céu aberto estão dando boas provas em todo o mundo hã já quase uma dezena 
de anos, a partir do incremento da morfologia matemática no campo mineiro (ver OBJECTIVO-INTRO- 
DUÇÃO deste trabalho). 

Na falta de experimentação real na aplicaçao do método agora descrito, limitamo-nos a apresentar 
o desenvolvimento do cálculo para um piso da Mina do Moinho (Aljustrel), a título de exemplo, 


PISO 200 (i = 200) 


O piso está dividido em 9 unidades de exploração que realmente se podem esquematizar do modo 
constante da figura 7. A 5.º unidade de exploração é uma unidade de franja composta por 3 câmaras. 


Poço de extracção 
Galeria de rolagem 


Camara da unidade de ex- 


ploração nº x 


Fig. 7 


Esquema simplificado do piso 200 da Mina do Moinho 


As características consideradas para cada unidae de exploração são as constantes do seguinte quadro: 


Unidade Reserva a Explorar | Distância média | Distância média | Teor de Teor de 


| de RB, | de rolagem (m) de carregamento enxofre cobre 
Exploração (ton)) d,, im) dj (% 8) (% Cu) 
1 32 700 295 25 + 25 = 50 46,91 0,62 
Pi 24 200 256 25 + 25 = 50 45,83 0,62 
3 38 500 2098 2) + 20 — 40 47,53 1,78 
4 24 800 342,5 25 + 25 = 50 47,07 1,80 
5 150 600 462,5 20 + 20 = 40 46,87 0,41 
6 58 600 315 33 + 33 = 66 45,74 0,83 
T 43 300 263 27 + 27 = 54 46,25 0,98 
8 T1 300 Sdo 26 + 26 = 02 47,51 1,54 
q 52 700 282 30 + 30 = 60 46,12 1,25 
Médias Ponderadas... ... ... D, = 347,8 L, = 49,8 46,68 0,93 


Reserva total a Explorar: R,,, = 5 R,, = 496 700 ton 


Deste quadro não consta o rendimento em peso, parâmetro fundamental normalmente para outro tipo 
de jazigo, porque se verifica que neste caso ele é igual para todas as unidades de exploração, o que se justi- 
fica facilmente por o minério ocorrer em «massas. 

Segundo o método exposto, as condições de exaustão do piso e de manutenção da frota constante 
implicam que o centro de massas do piso se mantenha constante. Então as possibilidades de combinação no 
piso resumem-se a combinar a unidade n.º 5 (única unidade de exploração com distância média de rola- 
gem superior à distância média ponderada de rolagem superior à distância média ponderada de rolagem de 
todo o piso) com todas as outras unidades de exploração. Estas combinações terão de se sujeitar, portanto, 
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à relação de 1 para 2,3 em tonelagem a extrair, visto que a unidade n.º 5 tem 150 600 toneladas exploráveis 
e o restante piso tem 346 100 toneladas exploráveis no seu conjunto. 
Segundo a simbologia usada, os valores reais que temos de considerar neste exemplo são: 


e,00 — 105 300 ton/ano = 0,545 ton/min, 
V, = 8km/h = 133 m/min, 
V,=K V,= 10km/h 
Vo 
K = = 1,25 
V, 
t, = 4 min 
Q' = 17 ton 
nv = 0,75 
V=0,5 m/seg = 30 m/min 
Q = 5 ton 
t = 6 min 
7 = 0,50 
q' = 1,9 ton 
V=34 m/seg = 2U4 m/min. 
T = 45 seg = 0,75 min, 
E IRA = 228 m 
Q =2q'=3,8 ton 
Trago — 0,80 


Considerou-se que um ano tinha 230 dias de trabalho e que cada dia de trabalho tinha apenas 
14 horas úteis efectivamente (o que é a realidade em dois turnos produtivos por dia). 
Com estes valores e por aplicação directa da fórmula (B) temos: 


N = 0,572 locomotivas 


200 


logo, uma locomotiva é suficiente para o trabalho de rolagem do piso. 
Os ciclos para cada unidade de exploração são: 


&, = 1,99 min 
Tc, = 7,19min 
e, = 7,49min 
c, = 8,63min 
z's = 10,26 min 
e, = 8,26min 
Eq — T,56 min 
cg => 8,67min 
zy= 7,82min 


A possibilidade de produção horária de uma composição para a unidade de exploração n.º 5 será de: 
P. — 99,4 ton/h 
mas como a taxa de desmonte daquela unidade é 
e, = 9,9 ton/h 
o número de composições por hora a trabalhar para a unidade n.º 5 é: 
N,., = 0,10 composições/hora 
Para as restantes unidades temos, com 
e, = 22,8 ton/h 
seguintes valores; 


= 0,18 composições/hora 
= 0,16 composições/hora 
= 0,16 composições/hora 
= 0,16 composições/hora 
= 0,18 composições /hora 
= 0,17 composições/hora 
= 0,19 composições/hora 
= 0,17 composições/hora 


”i “+ 
ta mer 


% 
a 


“j 
=, 
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—T 
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Note-se que a quase constância encontrada para todos estes valores de N, (exceptuando, logica- 
mente, o de Nes) vem provar a excelência do método. 
O número de vagons a encher também é facilmente calculado, Por exemplo, para a unidade n.º 5 


encher-se-ão: 


0,89 vagons/hora 


por aplicação directa da fórmula deduzida para o efeito. 
Cada composição será de 9 vagonas na unidade de exploração n.º 5, 


Temos, assim, totalmente resolvidos os problemas levantados pela 1.º equação do nosso sistema 
de equações que, representa o conjunto de relações estruturantes que definem as associações permitidas 
entre as unidades de exploração do piso considerado. No nosso caso, e dados os parâmetros constantes do 
quadro com as características de cada unidade de exploração do piso (os parâmetros são: distância de 
rolagem, distância de carregamento, teor de enxofre e teor de cobre), teremos um sistema de 4 equações 
a 4 incógnitas para resolver (5 = 4). Deste sistema já temos as combinações possíveis à luz da solução 
de uma das equações apenas, a equação referente ao parâmetro distância de rolagem (que foi a conside- 
rada «mais forte»), equação bem resolvida como se constata pelos valores de N, achados. Há agora que 
continuar a resolução do sistema de equações, tal como se expôs na parte final de 5. ESTRUTURAÇÃO 
DO MODELO, 


Como combinações possíveis para satisfação da 1." equação, recorde-se, tinhamos achado: 


1.* Unidade 2." Unidade 


ago 


e agora tomemos para 2." equação a referente aos teores de cobre (considerada mais fraca que a 1.*, mas 
mais forte que as restantes duas). O teor de cobre médio da unidade n.º 5 é inferior ao teor de cobre 
médio ponderado do piso, logo são de excluir as combinações 5-1, 5-2 e 5-6 e de considerar as combina- 
ções 5-3, 5-4, 5-7, 5-8 e 5-9, para as quais se obterão os seguintes valores médios ponderados de cobre: 


5-3 — 1,86 4 Cu (desvio de 46 %) 
5-4 — 1,44% Cu (desvio de 55%) 
5-7 — 0,81 9% Cu (desvio de 13%) 
5-8 — 1,06 & Cu (desvio de 14%) 
5-9 — 0,99% Cu (desvio de 6%) 


sendo as ponderações feitas em função das taxas de desmonte de cada unidade de exploração (calculadas 
na resolução da 1.º equação do sistema) e os desvios referem-se à média ponderada do teor de cobre médio 
do piso. Todos estes desvios são no sentido positivo (percentagens de cobre superiores à média ponderada 
do piso), exceptuando o referente à combinação 5-7. As combinações excluídas são todas excluídas por 
combinarem duas unidades de teores médios de cobre inferiores ao teor médio ponderado de cobre do 
piso, pelo que todas essas combinações também terão teor médio ponderado inferior ao teor médio do piso. 
Assim, é lógico que não se deva considerar todas as combinações achadas possíveis até agora e se deva 
fixar um y para desvio máximo admissivel positivo (teor médio ponderado da combinação superior ao 
teor médio ponderado do piso), Se este y for 15% teremos de excluir as combinações 5-3 e 5-4 e con- 
siderar apenas as combinações 5-7, 5-8 e 5-9, 

Posto isto, passa-se para a resolução da 3.º equação do nosso sistema. É a equação referente ao 
teor de enxofre e o seu tratamento é em tudo semelhante ao da equação anterior. Só teremos como 
possíveis as combinações 5-7 e 5-9, para as quais se obterão os seguintes valores médios ponderados de 
enxofre: 


o-1 — 46,43% S (desvio de 1%) 
9 — 46-35 % 5 (desvio de 1%) 
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Com tão pequenos desvios em relação à média ponderada do teor de enxofre médio do piso, acei- 
tam-se as duas combinações, E numa 2.º análise pode-se também aceitar a combinação 5-8 excluída, por 
tê-lo sido por ambas as suas unidades de exploração terem teores médios de enxofre superiores ao teor 
de enxofre médio ponderado do piso. 

Passar-se-ia agora para a 4.º e última equação do nosso sistema e que é referente ao parâmetro dis- 
tância de carregamento. O seu tratamento é idêntico, pelo que nos dispensamos de o descrever, excepto no 
que diz respeito ao cálculo das necessidades de equipamento de carregamento. 

Para obtenção de tudo o que se considerou na exposição do método descrito, basta aplicar a fórmula 
(A) directamente e teremos: 


Doo — 1,670 pás 


logo são suficientes 2 pás auto-carregadoras no piso, distribuindo-se uma para cada unidade de exploração 
em actividade. A produção horária deste equipamento de carregamento é igual à produção horária já cal- 
culada para a frota de rolagem. 

Chegámos, assim, ao fim da resolução do sistema de equações, apenas com duas combinações a satis- 
fazê-lo totalmente, nas condições que se expõem no quadro: 


Tempo de es- RESERVA DESMONTADA 


Taxa de des- Reserva brméáto dá Reserva 
monte | inicial sia dad mind sobrante 
(ton/h) (ton) dades (dias) | (ton) Ton 
8,8 150 600 — 131 799 
62 101 
22,8 43 300 136 — 
9,9 131 799 =— 108 916 
T5 583 


22,8 2 700 165 a 


Este quadro pode e deve entender-se como o princípio de uma das sequências possíveis de explo- 
ração. Mas ele afasta muitas unidades de exploração de forma talvez demasiado rígida, como o é o facto 
de não se ter admitido desvio algum possível para o parâmetro distância de rolagem, que era o pará- 
metro da equação mais forte do sistema (a 1.º equação a ser resolvida). Para se ter uma noção da influén- 


cia desta fixação de desvio nulo, basta reparar que se admitirmos que D, se pode deslocar 5% para 
qualquer lado teremos um y = 17,4 metros ou seja D, poderá estar compreendido entre. 


D, +y= 347,8 + 174 = 3652m 
e D, — y = 347,8 — 17,4 = 830,4m 


e ao fazermos isto vamos colocar dentro deste inter valo as unidades de exploração n.º 4 e 8, as quais 
passarão a poder combinar-se com as outras unidades de exploração, para além de se poderem combinar 
entre si e de poderem continuar a combinar-se com a unidade n.º 5. Assim teriamos como possíveis mais 
estas combinações: 


4-1, 4-2, 4-3, 4-6, 4-7, 4-9, 4-8, 8-1, 8-2, 8-3, 8-6, 8-7 e 8-9 
para além das já consideradas: 
dl, 5-2, 5-3, 5-4, 5-6, 5-7, 5-8 e 5-9 


o que perfaz um total de 21 combinações, contra o anterior total de 8 combinações, só por se ter admi- 
tido o desvio de 5% num parâmetro que estava totalmente fixo. 

Ter-se-ia agora de resolver as restantes 3 equações do sistema considerando estas 21 soluções da 
1.º equação e, no final, completar o quadro da sequência possível de exploração. 

Para se estender o exemplo até à totalidade do campo estudado neste trabalho resta averiguar o que 
se passa no poço de extracção, Por aplicação directa das respectivas fórmulas temos: 


Zponn Em 2,94 min. 


Prsgo — 1,08 ton/min. = 61,8 ton/hora 
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e o poço tem apenas necessidade de trabalhar meio-dia para este piso em cada dia de trabalho (a taxa 
de desmonte era de 32,7 ton/hora). Não nos esqueçamos que esteamos dentro do caso habitual (Q' > Q), 
que chegarão à boca do poço 2 composições por hora (6 = 2), pelo que teremos no poço 18 vagonas por 
hora e a existência total de vagonas no piso terá de ser; 


I8 x (2+ 0) = 72 vagonas 


9. NOTAS FINAIS 


O exemplo apresentado não permite avaliar a total generalidade do método exposto, dado o pequeno 
número de unidades no piso estudado. Mas teve de ser assim por não dispormos, no local de trabalho, 
de capacidade de cálculo computorizado que permitisse a introdução de mais unidades de exploração e de 
outros pisos. Desejamos notar de novo que qualquer ordenador está apto a proceder a toda esta sequência 
de cálculo, depois de devidamente programado, dando em «out-put» uma sequência de exploração possível, 
a frequentar, depois, arbitrariamente, e os dados adicionais desejados. 

O piso dado como exemplo permite desde já afirmar que a sua actual situação é de desiquilíbrio, 
pois se ele fosse explorado segundo o método exposto não estaria hoje na situação em que ele se encontra 
e que pode ser caracterizada, por exemplo, pelo facto da sua unidade de exploração n.º 5 ter duas das suas 
câmaras em 8.º —«corte ou planta» e a outra câmara em «6.º corte» (recorde-se a importância funda- 
mental desta unidade de exploração no piso). Este desiquilíbrio ainda pode ser corrigido, se se proceder 
tal como se fez ao longo deste trabalho, pois os dados considerados são os reais nesta altura. 

Os autores agradecem à Administração de Pirites Alentejanas S. A. R. L. a autorização concedida 
para a apresentação do exemplo. Igualmente agradecem idêntica autorização da Administração de Mines 
et Industries S. A. embora esse exemplo não esteja tratado aqui. 


EM TEMPO: 


I. O presente trabalho, ao ser apresentado no Seminário de Geomatemática, foi precedido de algu- 
mas palavras introdutórias que focaram essencialmente o objectivo da «parametrização», por blocos tecno 
lógicos, em termos de projecto, e as vantagens de critérios de rejeição se basearem nos parâmetros «teor 
de concentrado» e «rendimento em peso». Igualmente foi feita referência à vantagem da estacionarização 
preferencial daqueles dois parâmetros, cuja variância é menor do que a referente aos parâmetros do tal- 
-qual, tal como foi devidamente demonstrado pelo Sr. Prof. Eng." Quintino Rogado no seu trabalho apre- 
sentado ao 13% APCOM, 

Na sequência destas observações foram tecidas considerações sobre as vantagens da estacionarização 
de parâmetros e das consequências de fixação de um «teor de cortes, nomeadamente para: 


-—— saber como evolui o teor médio global ao longo da vida da Mina; 

— saber se há necessidade de alterar frequentemente o teor de corte ao longo da vida da Mina; 
— garantir o teor do concentrado vendável; 

— definir as condições de «blendagem» na Mina; 

— garantir nas reservas sobrantes, em cada instante, as características iniciais do jazigo, 


— actuar na definição das caracteristicas contratuais de venda dos concentrados e assumir os res- 
pectivos riscos de penalidades conscientemente; 


etc., etc. 


II. No final, após a exposição do trabalho, foi apresentado um outro exemplo, referente ao piso 265 
da Mina do Moinho (Aljustrel), para completar o exemplo constante do «pré-print» distribuído aos partici- 
pantes no Seminário em causa uma vez que este piso tem substanciais diferenças em relação ao outro e, 
também, porque se entendeu conveniente fazer um estudo de interligação dos dois pisos, devido a pro- 
blemas actuais que se sentem naquela Mina referentes à utilização do poço de extracção, 
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Dado o pormenor a que se desceu no «pré-print» na apresentação do exemplo, limitamo-nos agora 
à apresentação dos dados iniciais e das conclusões referentes a este novo exemplo. Assim: 


— Esquema do piso: 


Para o poço 
de extracção 


Fig. 8 


Esquema simplificado do piso 265 da Mina do Moinho (simbologia idnêtica ao do esquema da fig. 7) 


Note-se: 


a) a influência da falha no traçado das galerias de rolagem e na existência de 3 câmaras nas 
unidades de exploração n.º 5 e 6; e 


b) a existência de 3 unidades de franja (duas com 3 câmaras cada e uma de uma só câmara). 


— Características de cada unidade de exploração: 


Unidade Reserva a explorar Distância média Distância média Teor de | Teor de 
de BR. de rolagem de carregamento enxofre cobre 
exploração (ton) d,. im) hi. (m) xs c& Cu 
1 323 700 280 19+19=38 46,86 0,81 
& 310 300 | 235 26+26=52 45,95 0,70 
3 130 400 265 25+25=50 46,09 1,37 
1 225 700 193 15+15=30 46,31 1,58 
5 467 200 225 20 + 20= 40 | 48,08 1,21 
6 359 600 367 22+22=44 | 47,96 1,05 
RES DR E O E as Es e. 
Médias ponderadas... ... ... L =42,0 47,11 | 1,01 | 


Reserva total a explorar: R,,,. == R, = 1816 900 ton. 


— Repartição da reserva em função de DB: 


R(+) ... 8137 00ton ... 44,8 % 
R(—) ... 1003200ton ... 55,2 % 


— Valores a usar são todos os do piso 200, excepto: 


Eis — 134 000 ton/ano = (,694 ton/min 
Oras (+) = 0,811 ton/min = 18,86 ton/h 
e (—) = 0,383 ton/min = 22,98 ton/h 


20% 


total -s 0,694 ton/min = 41,64 ton/h 
H,,; — 288m 
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— Resultados obtidos: 


Na; — 0,412 locomotivas no piso 


Z, = 7,79 minutos 
2, = 7,18 minutos 
Z', = 7,59 minutos 
z', — 6,61 minutos 
z'; = 7,05 minutos 
z'; = 8,97 minutos 
P. = 130,9 ton/h 


P, = 1421 ton/h 

P, = 134,4 ton/h 

P e 154,3 ton/h 

P, = 144,7 ton/h 

P. = 113,7 ton/h 

Ne, = 0,14 composições/hora 

N.., = 0,16 composições /hora 

N x = 0,14 composições/hora note-se esta quase constância obtida 
N,, = 0,15 composições/hora | (é um bom resultado) 
N.., = 0,16 composições/hora 

N,. = 0,16 composições /hora 


1,35 vagonas/hora em cada unidade de exploração 


19 vagonas/dia em cada unidade de exploração 
2 composições/dia em cada unidade de exploração (um comboio de 10 vagonas e outro de 9 vagonas) 


Do, — 205 =2 pás no piso 


z = 3,57 minutos 
P265 
= 0,85 ton/minuto = 51 ton/h 


Prego 


e, = 41,684 ton/h, o que equivale a cerca de 80% de um dia de trabalho do poço de extracção 


Q>| 
g = 


19 X (2+ 0) = 76 vagonas no piso 


— Início de uma possível sequência de explora ção 


Tempo de esgo- Reserva desmontada 


Taxa de des- Reserva Reserva - 
Combinatória monte | inicial 6 a a sobrante pómalar 
(ton/h) (ton) (ton) 
(dias) 
1 18,66 323 T00 1239 — 
5 22,98 467 200 — 68 560 
6 18,66 S09 600 ema 107 633 
2 22,98 310 300 865 — 


Note-se que se admitiram os mesmos desvios considerados para o piso 200, ou seja D, rígido, 15 % 
para o cobre e 1% para o enxofre. 


— Observações pertinentes: 


Não repugna o sobrecarregar uma composição de vagonas com uma vagona mais do que o habitual 
(10 em vez de 9 vagonas), mas há que notar bem que n,,. veio a dar 2,05 pás como necessárias para Oo piso 
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em estudo, Perante este número põe-se a questão: teremos de usar 2 ou 3 pás no piso? A resposta lógica 
seria 2 pás, para o que bastaria melhorar um pouco o rendimento de cada uma delas, mas a prática mos- 
tra que, com tal número, se opta sempre por colocar no piso uma 3.º pá auto-carregzadora, a qual ficará 
em regime de sub-utilização acentuada, embora permita uma melhor manutenção das 3 unidades em tra- 
balho no piso. 


— Relacionamento entre os pisos 200 e 265: 


Piso 200 Piso 265 
Csa0o — 32,7 ton/h Coogs — 41,64 ton/h 
Prov = 61,8 ton/h Prog; = 51,0 ton/h 
tempo de trabalho do poco = 50 % do dia idem = 80 % do dia 


O dia de trabalho sempre considerado foi de 2 turnos, com um total de 14 horas diárias produtivas, 
mas o já apontado leva a uma utiliação de 50 + 80 = 130 % daquelas horas num total máximo possivel de 
100 %, o que é manifestamente impossível. Uma solução é usar 50% para o piso 200 (a totalidade neces- 
sária para o piso) e outros 50% para o piso 265, ficando 30 % de 14 horas para ser extraído do piso 265 
em 3.º turno de trabalho do poço, o que corresponde a 4,2 horas de trabalho (teríamos sem trabalho apenas 
7— 4,2 =2,8 horas úteis para possível manutenção). Esta solução, indesejável, tem ainda o inconveniente 
de exigir um volante de vagonas no piso, subindo substancialmente o número de vagonas necessárias no 
piso. Para que o poço trabalhe apenas 2 turnos por dia põôem-se três alternativas: 


— melhorar o rendimento r, (qualquer melhoria possível deve ser insuficiente para o desejado); 
— ampliar a capacidade de extracção do poco (o que exige um investimento adicional), ou 


— fazer transbordo intermédio algures, para outro poço de extracção, de dada percentagem do miné- 
rio a extrair do piso 265. 


Verifica-se, assim, que nesta Mina, contrariamente ao suposto por nós em 7, o poço de extracção é 
«gargalos da produção para os níveis inferiores do piso 200. Neste caso concreto foi escolhida a 3.º alter- 
nativa das anteriormente apontadas, cuja extracção se esquematiza, no que se refere ao excesso de capa- 
cidade do poço que se tem apontado, do seguinte modo: 


POÇO POÇO 


Piso 150 
Piso 200 
Piso 265 


Fig. 9 


Esquema em vigor na Mina do Moinho (elevação do piso 265 até ao piso 150, aonde se usa uma 

chaminé-silo colectora para fazer chegar o minério a uma galeria de rolagem, do piso 200, com 

ligação a um outro poço de extracção que já era pré-existente, pelo que apenas se investiu na 
chaminé-silo colectora) 
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tContinuação da pdg.º 260) 


nuclear húngara, projectada em 
1966, localiza-se em Paks e será 
equipada com quatro reactores 
com uma potência total de 1760 
MW (tipo de reactores PVR), de- 
vendo o primeiro entrar em funcio- 
namemnto em Dezembro de 1980 e 
os restantes três em 1985. O equi- 
pamento para esta central, cujo 
custo total está calculado em 26 
mil milhões de forints, será forne- 
cido pela URSS e pela Checoslo- 
váquia, trabalhando no empreendi- 
mento técnicos soviéticos e pola- 
cos. 

Polônia — Em 1974 a Polónia e a 
URSS assinaram um acordo com 
vista à construção da primeira cen- 
tral nuclear na Polónia, seguido, 
em 1976, pela conclusão de um 
acordo de cooperação relativo à 
produção de instalações para cen- 
trais nucleares, A central será lo- 
calizada em Zarnowiec, no litoral 
de Gdansk. 


O início dos trabalhos deverá 
ocorrer em 1979 e a entrada em 
funcionamento da central em 1984, 
Na primeira fase, a Zarnowiec 
será dotada com dois reactores 
VVER-440 MW, aos quais se acres- 
centará, mais tarde, um terceiro, 
com a potência de 1000 MW. Em 
1988 a potência da central atômica 
polaca elevar-se-ã a 1880 MW, e a 
capacidade nuclear do país será se 
4880 MW em 1990 e de 23 000 MW 
no ano 2000, 

República Democrática Alemãd-— 
A RDA possui duas centrais nu- 
cleares em funcionamento, a de 
Rheinsberg, ligada à rede desde 
1966, e a de «Bruno Leuschner 
Nord I>» em Lubmin (na baía de 
Greifswald) desde 1973. Esta últi- 
ma atingiu em Maio deste ano 
uma potência instalada total de 
1520 MW com a entrada em ser- 
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viço, com o auxilio soviético, do 
terceiro reactor de 440 MW. Pre- 
sentemente, a RDA procede à am- 
pliação da central de «Bruno Leus- 
chner Nord II», com vista a atingir 
uma capacidade instalada de 3520 
MW. 


Checoslováquia — Desde ISTO, 
data da assinatura do acordo inter- 
governamental entre a Checoslová- 
quia e a URSS sobre construção 
de centrais nucleares, a Checoslo- 
váquia vem a realizar, com a co- 
laboração da União Soviética, o 
seu programa de construção de 
centrais atómicas. 

Em Jaslovske Bohunice, na Es- 
lováquia, encontra-se já em funcio- 
namento desde 1973 uma central 
experimental AI (urânio natural, 
água pesada, gás) de 143 MW, es- 
tando a ser construída, no mesmo 
local, uma segunda central VI do- 
tada de dois reactores VVER-d40, 
com uma capacidade individual de 
440 MW. A entrada em funciona- 
mento dos dois reactores está pre- 
vista, após vários atrasos, para 
1978 e 1979, respectivamente. Em 
1976, foi dado início aos trabalhos 
de edificação de uma terceira cen- 
tral V2 com características seme- 
lhantes às da VI, em Dukovany 
(Morávia do Sul). O primeiro 
reactor começará a funcionar em 
1982, ou seja, com dois anos de 
atraso, e o segundo em 1984. A 
sua construção terá a assistência 
técnica soviética. 


Em 1985, a parcela representada 
pela energia nuclear na produção 
energética total atingirá 17 %, em 
1990 cerca de 40 % e no ano 2000 
cerca de 70%. A concretizarem-se 
as previsões, a Checoslováquia de- 
verá possuir uma capacidade nu- 
clear instalada de 800 MW em 
1980, de 3500 MW em 1984, 10 500 
MW em 1990 e de 30-35 000 MW 
no ano 2000. 


A primeira central (experimen- 
tal) da Checoslováquia foi reali- 
zada pela URSS. Relativamente à 
segunda, em construção, a partici- 
pação soviética limitou-se ao for- 
necimento de determinados ele- 
mentos do circuito primário. De 
futuro a indústria checa devera 
construir na integra as centrais a 
instalar no seu território e forne- 
cer igualmente determinados equi- 
pamentos à URSS e demais países 
socialistas. O país especializar-se- 
-à, nomeadamente, na produção de 
reactores do tipo VVER a água 
pressurizada, prevendo que a sua 
produção anual, parcialmente des- 
tinada à exportação para os pai- 
ses do Comecon, atinja as cinco 
unidades antes de 1980. Este tipo 
de reactores será produzido até 
1990; por volta de 1985, a sua 
capacidade será elevada para 1000 
MW. A produção do primeiro reac- 
tor a neutrões rápidos e arrefeci- 
mento por meio de sódio foi já 
iniciada, 

Por outro lado, encara-se a pos- 
sibilidade de estabelecimento de 
novo acordo entre a Checoslová- 
quia e a URSS, mediante o qual 
este pais forneceria ao primeiro, 
electricidade contra a entrega do 
equipamento checo para a sua 
central de Khmel'nickij, bem como 
de diversos bens de consumo, 


Cuba — Em 1977, Cuba iniciou, 
com o auxílio da União Soviética, 
a construção da primeira central 
nuclear, que deverá estar conclui- 
da em 1980, na região de Cienfue- 
gos, 300 km ao sul de Havana. A 
URSS fornece os equipamentos ne- 
cessários à Central, cujo primeiro 
bloco terá uma capacidade de 400 
MW, estando a capacidade total 
do complexo calculada em 1700 
MW. 


[Extraido do «Boletim Mensal BFNs», 
(12) Dezembro 1978]. 
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Process Analysis and Simulation in Extractive 
Metallurgy by Stochastic and Deterministic Methods 


ALEX OBLAD 


Neste artigo apresentam-se métodos correntemente usados na construção e análise de modelos empi- 
ricos pela via estocástica clássica e também recorrendo a modelos determinísticos recentemente criados. Nesta 
última óptica e usando o modelo de balanço generalizado em sistemas de partículas, tratam-se aplicações 
à moagem discontinua/continua, à flutuação e à lixiviação. 


ABSTRACT 


This paper deals with current methods used in construction and analysis of empirical models by lhe 
classical stochastic approach and with recent developments in determiníistic models. These latter include the 
applications of a generalized population balance model for particulate systems to the operations of batch/ 
“continuous grinding, flotation, and concentrate leaching. 


INTRODUCTION 


The optimization of mining methods often is influenced by the processes that will be used to treat 
the ore. The behavior of these processes with respect to operating conditions and material feed characte- 
ristics is of prime interest, and considerable effort has been expended to define this behavior, not necessarily 
for applications to mining methods, but in general to gain a detailed knowledge of how to control the process 
products as a function of the presumably varying mineralogical characteristics of the feed. 


In accomplishing this, simulation of unit operations and more recently entire plants has received the 
most direct attention. Process simulation has been approached from two contrasting methodologies. The pur- 
pose of both, however is to develop mathematical models for a given process. 


The first approach, empirical modeling, samples the process response at various points in the space of 
the variables considered to control the response, and fits a surface of predetermined form to the sampled 
data. This approach has little or no theoretical basis, and the parameters, i.e. constants, are not physically 
meaningful. The model is entirely dependent on experimental measurement, and is therefore composed par- 
tly of real process variation and partly of error. Extreme rigor must be employed to assure that the result 
is a model essentially of the true process response and not that of the error. 


Deterministic models are considered to comprise two types: phenomenological and theoretical. In the 
first, the equations for the process response have a theoretical basis but include parameters whose values 
can not be determined without experimental data. These parameters in general have physical significance. 
The theoretical model completely specifies the process response with parameters which can be calculated 
from theoretical considerations. Phenomenological models have achieved some success in industrial applica- 
tions. Unfortunatelyv, theoretical models usually require too detailed information to be of practical use, at 
least in extractive metallurgy. 


The use of both stochastic and deterministic methods is now receiving almost universal application 


in the minerals industry. This paper will discuss some of the more significant developments in certain areas 
of mineral processing and hydrometallurgy. | | 
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STOCHASTIC METHODS 


In the context of process analysis, stochastic modeling involves the application of classical statisti- 
cal inference to estimate model parameters, to determine intervals in which the true values of such para- 
meters lie, to determine the probable error of both data and predictions; and to provide a basis for testing 
of hypotheses about the model parameters and responses. 


In general, these models are empirical. This precludes the evaluation of parameters in terms of physi- 
cal reality; that is, the model parameters have no physical meaning. The experimental data are fit to a 
response surface in the space of the independent variables, the model having little or no theoretical basis. 
It is therefore extremely important to know the relationships between the error associated with the expe 
rimental measurement of the phenomenon being studied. This, of course, is true for all types of models, but 
especially for empirical models. 


EMPIRICAL MODELS 


These models are classified into two main types: linear and non-linear, Linear models are those which 
are linear in the parameters, or are transformable to linear. For exemple the equation: 


Y=1-exp (=Kx”) 
is transformable linear: 
Z=atbW 


where Z— ln (—ln (1-y)) 


W= In X 
a=ztm K 
b=N 


To estimate the parameters of any model the leastsquares criterion is usually invoked. If we pos- 
tulate that Jf(K, X,... X: Bo 8, B,) in which the X are the independent variables and the p, are 
parameters, then the measurement of the deviation between the predicted and experimental response, called 
the objective function, will be defined in the context of the least-squares criterion as follows: 


O = ba bió (o hi E (XX AM 5 Bos B By) j 1 
where 4 = objective function 


Y, — experimental response for experimental level j 
; = predicted response for experimental level j 


w, = weighting factor for experimental level j 


The objective function q will be minimized by choosing the appropriate values of the parameters 
Bo B, o B, If f is linear in the parameters, then an analytical solution is possible. Otherwise, successive 
approximations must be used based on some type of implicit linearization, Only linear functions will be di 
scussed here, It will also be assumed that the independent variables are non-stochastic, i.e. that they can be 
specified with negligible error, 


Let the true response of a linear system be represented by: 


n=B+ 8 Ny) + BA — Ng) + cce BAN, —N,) é 
Each experiment will correspond to a different level of the independent variables except when there 
are replicate experiments made at a constant experimental level to estimate experimental error. 
It is possible to define a design matrix X whose elements *|j — x, correspond to the value of variable i 
at level j. The vector of experimental responses Y has elements y, which is the average response at level j; 
and let B be the vector of elements 8, Im this formulation 


do 2 (PD; *y) I 2 P; 3 


where Pp = number of replicates at level j. 
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Therefore, g =— > Wj ( Y; — ny) 4 
j=i 
or in matrix notation: 


[ = pe W [E =] 5 
where 
E = 
W 1j ! , =X Ê 
o Ysk) = 
To find the vector 8 giving the minimum value of g, it is necessary to solve the equation: à 9/0 p = o. 
0 QNT f F W 6 
E E (Sep 
XWY=Xw 8 7 


E [XUW XT exists, then the best estimate of / is 

&=[X"N XP XWY 8 
Therefore, the best estimate of q is: 

P=X8=X[X'WX|PXW yY 9 


The variance and covariance of the parameters can be calculated from: 
Cr=[XWX7]'S 10 


where C is the variance-covariance matrix and q? can be estimated from Ses 


* o | dr] Pa e = F ua 
C,7=VAR(B,,) C,7=COVAR(É,| 8.) 


For the predictions: 


E e oe ”, 
where d7=VAR(?,), “d =COVAR(Í, 7) 


To test the hypothesis that the model describes the system within experimental error, compute 


N Pi ii 
| EE 


— aa 


2 pj— N 
N o a ol 
2 ds Yi) 
 N-00 13 


where Q = number of model parameters, N = number of levels. 


“Bl 


12 


If Sº/B < F. (N-Q, = P,—N), then the model is judged to be correct within the limits of expe- 
E 
rimental error. If the model is correct, then an improved estimate of 5? would be: 


E ga ris TO E A qua o VE 
oo 2 vily=- A) + 2 23 (Yn- 44) dá 
2 PO 


To determine joint confidence intervals for the model parameters, the solution to: 


b ! “+ 


[B-BI XWxX( B-8 )<Qo Ra 

defines the (1-w) 100 % confidence limits in parameter space, With the appropriate experimental planning, 
for example the 2-level factorial composite designs, empirical models can be developed from experimental 
data, and probable variations of the parameters and predictions can be assessed relative to experimental 
error. This gives a quantitative measure of whether the model is a significantly good representation of the 
data, and a basis for choosing between of competeing functional forms. 
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DETERMINISTIC METHODS 


Extractive metallurgy deals in general with particulate processes involving necessarily a countable 
number of entities. These possess a distribution of physical and chemical properties, the most common of 
which is the particle size distribution. 

Distributions are expressed mathematically by either the cumulative distribution function F (Z) or 
the density function f (z). F (£) gives the number (i=1), área (i=2) or mass (i=3) fraction of the particles 
with property z < é. The differential quantity f (£) dz gives the corresponding fraction of the particles with 
property z in the range £ to : + dz. 

Mathematical modeling of particulate systems from the deterministic viewpoint attempts to develop 
general theoretical approaches to be used in determining the detailed kinetic behavior of the appropriate 
property distributions. To accomplish this it is necessary to consider, besides mass, momentum and energy 
balances, a population balance which conserves and identifies particulate distributions. In the next section 
the mathematical structure of such a balance is developed. 


POPULATION BALANCE MODEL 


To construct the population balance model we must define the following properties: 


Let: 
n(xX, y. Z,t) = Total number of partícles per unit geometrical volume at position x,y, Z and time t. 
f(£; x,y,Zz, t) dZ = the fraction of particles at position x,y, z and time t with property Z in the range 
é to E + dZ. 


v=n(x,y,21t) f,(£ x,y,2,t) dZ dz dy dz = the population density: lhe total number of particles 
in the region of geometrical space |[x, x+dx], [y, y+ 
+dyv], [z, z+dz] and in the region of «property» space 
[£, t+dZ] at time t. 


u=da/dt = the property kinetics. 

B(£;x,y,2z,t) dZ — the number of particles at position x, y, z and time t created per unit geometrical 
volume per unit time with property in the range £ to £ + dZ. 

D(£;x,y,2z,t) dZ = the number of particles at position x,y,z and time t destroyed per unit geometri- 
cal volume per unit time with property Z in the range £ to E + dZ. 


The quantities B and D are called the birth and death functions, and account for changes due to ins- 
tantaneous appearance or disappearance of particles, as oppossed to continuous changes. A familiar example 
would be the distribution of particles in grinding due to breakage in which the resulting daughter fragments 
are suddenly created. 


A population balance requires that [1]: 


Oy d 
—— + N.(vy) + (uy) +D-B=O 
at dZ 
dy 
where —— = change in population density per unit time, 
dt 


V .(vy) = change (divergence) in population density due to motion in geometrical space (v is geo- 


metrical velocity). 
| 


(uy) = change in population density due to the continuous kineties of property, 
Ef 


D—B = net change in population density due to instantaneous creation and destruction of particles. 
In most applications it is not convenient nor even possible to specify the detailed spatial dependence 


of y; thus, volume-averaged population densities must be used. The resulting equation is called the macros- 
copic population balance model: 
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Ito 7) e: PN ga 

V,dt + v (qa da) np )+D=B= Z 
where 

e A 4 o e RE] | o | o 

b = Er | vdV, D= = DdV, B= V BdV; QO,, OQ, = volumetric flow rates 


at reference points 1, 2; V = volume considered. 

It is assumed that the knetics of property Z is position independent, so that u = u(t). 

Equations 1 and 2 present at least in principle a detailed and precise means of describing quantita- 
tively the behavior of particulate systems in terms of property distributions. This is especially appropriate in 
mineral processing operations because products are generated whose properties are commonly measured with 
respect to distributions of particle size, grade, density, and other variables. In the most general formulation 
the population balance model is not analytically solvable, but for several important process applications it 
has been possible to simplify the equations to the point where meaningful solutions can be obtained. The pre- 
dictive capability of these models is superior to empirical or less rigorously derived equations. The remain- 
der of this paper will discuss applications of the model to the operations of grinding, flotation and leaching. 


GRINDING 


Historically, three approaches have been used in the analysis of tumbling mill grinding. These are 
Energy Size-Reduction relationships, Matrix Models, and Population Balance Models. 

Energy size-reduction relationships ignore all details of feed and product size distributions and are 
empirical in nature. These relationships have been oriented towards evaluating the amount of energy needed 
to achieve a certain value of a single Kinetic parameter chosen somewhat arbitrarily to represent the pro- 
duct size distribution. Since tumbling mills draw power at a constant rate, energy may be replaced by 
time to transform the relationship into a single-parameter Kinetic expression. 

The majority of these relationships can be written in the general form: 

dEv 
dX 


=— KX* ú 


where E = Energy input per unit volume applied to the particulate system, 
* = A characteristic particle size of the system, 
The kKinetic relationship is obtained by substituting Ev = Pt/V, 
where P = power draft, V, — volume of solids, 
dX [E PR 
dr o K Vs 


For evaluating energy needed to achieve a certain value of X, integration of equation 3 yields: 

KR [os IN | 

—— | Ao o — A. + Nasd 

Ev d N-] l e | Lá 

K In (X/Xo) a E 

where x. = value of X in feed or when t= o depending on whether the grinding is continuous or batch. 
For different values of N we have the famous «Laws» of comminution. 


Kick's Law (1885): Equivalent amounts of energy applied to the system should result in equivalent 
changes in the size of the material being ground. This implies equal sized fragments and breakage resistence 
independent of particle size. Equation 5 becomes for N=1: 

E, =K, In (X/X,) 6 


Rittinger's Law (1867): energy consumed is proporcional to the new surface produced per unit mass. 
Equation 5 becomes for N = 2: 


| l 
Ev = ks X ese | 


Bond's Law (1952): energy consumed is proportional to the total length of new cracks formed. Equation 
5 becomes for N = 3/2: T 


a I 
Ev = K q— | — 8 
: o VA VX, ) 


where X, X are the 80 % passing size of product and feed, respectively. 
Bond's Law has been applied with considerable success in the dimensioning of industrial mills, This 
is due to the fact that equation 8 is satistied if the batch grinding product size distribution in the region 
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of application has a slope of 1/2 in the appropriately linearized transformations of the Gaudin-Schuhmann 


or Rosin-Rammler distributions. Fine grinding often has batch size distributions fitting these conditions. 
To approximate continuous grinding, multi-batch tests can be made until the mill product ceases to change. 
This state of equilibrium is equivalent to a continuous plug-flow mill operated in closed circuit. Since par- 
ticle residence time distributions in mills approximate perfect mixers in series [2], this assumption is an espe- 
cially bad one. Nevertheless, Bond's correlation is the most common method to perform mill calculations. 


The second approach, that of Matrix Models, proposes a matrix D which transforms the feed size 
distibution to that of the product for each breakage event as follows: 


ki 
J 


where W,, = Mass fraction of product in size fraction i 
Wp, = Mass fraction of feed in size fraction j 


J 
D,, = Mass fraction of feed size j appearing in product size i 


The matrix D is divided into two parts. Let 8, = the mass fraction of material in the ith size class 


which is broken in one breakage event, and B; = the mass fraction of the breakage product from the ith 
size class which enters the jth size class during a breakage event. 
A material balance after one breakage event yields: 


i 
W ai E (1-Sf) Wp + >, Bi Sj We = 
d 
In matrix notation equation 10 becomes: 
Wp=[1—=S+B'S]W,=DWP u 
where 
“40 0 
UU |) 
VOL) nua sm 
[= PO Dis sb RR QAEDA, TARLVA 
00001 O 
O Q00 01 Ú 
O OO | OO La 
se P 
sg 
SE == Ss. nxn diagonal matrix 
U 
Er 
30 nt 
Do, Bag |) 
Bê= Re nxn lower triangular matrix 
B . B; 


E o DR RS PG E 
For v identical stages of breakage: 
Wyi=D" W, 12 


To arrive at an expression for the kinetics of grinding, the cycle factor v must be related to the real 


time variable. Once this is done the cycle selection functions sº are estimated from experimental data; 
and in practice it is assumed that 


N - 
D tu =| 1=esp(aijas) |/[ 4 esp | 13 


|) 
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a modified form of the Rosin-Rammler distribution. The variable d represents particle size. 

Matrix models have been found to be useful in simulating single-pass devices in which the contact 
time is on the order of the time needed for particle fracture. For dynamic systems these models are inap- 
propriate. 

In contrast to Energy-Size Relationships or Matrix Models, Population Balance Models have the follo- 
wing advantage: 


1) They represent the behavior of the entiresize distribution in contrast to Energy-Size Reduction 
correlations which characterize the size distribution with one parameter, 

2) “The time variable is explicitly accounted for,in contrast to Matrix Models which must correlate 
breakage cycles with time; 

3) Process behavior is presented in terms of physically meaningful quantities associated with the 
actual mechanisms causing size reduction; 

4) The averaging process inherent in the macroscopic population balance model effectively allows 
the behavior of the grinding process to be specified without a detailed knowledge of the mecha- 
nisms involved. In this sense, part of our ignorance about the spatial details of grinding is 
suppressed without severly prejudicing the model. 


For batch grinding equation 2 can be arranged for solution from the following considerations: 


Q,.=Q,-0 14.a 
Z=d 14.b 
d(d)/dt=o 14.c 


Equation 14.a signifies that the system is batch; equation 14.b, that the property of interest is the 
particle size; equation 14.c, that all continuous changes (for example, attrition) are negligible. 
The macroscopic population balance model reduces to: 


à sor no a 
| n(do) |=B-D 15 


where Dd(d) is the rate of breakage (number/volume-time) of particles having linear dimension d to d+d(d); 
Bd(d) is the rate at which fragments, produced from breakage of sizes greater than d, enter the size inter- 
val d to d+d(d). 

Define S(d, t), the Selection Function such that: 

S(d,t) n f(d,t) d(d) dt = number of particles of size d to d4-d(d) broken in time interval t to t+dt; 
and b (d, d'), the number Breakage Function, such that: 
b (d, d') d(d) = number fraction of the fragments produced from breakage of particles of size d' 


1 


which enters the size interval d to d+d(d). 

The rate at which particles of size d' are broken is: 

S(d',t) = af, (d',t) d(d'), and the contribution to the interval d to d+d(d) resulting from the brea- 
kage of particles of size d' would be: 


The total contribution from all sizes d' > d is: 


“dd max 
B d(d) = | b, (dd) S(d't) q f, (d'st) d(d”) d(d) 16 
d 


Equation 15 thus becomes: 


“A a a 
e | nf(d,t) | = — S(dyt) n fi(dist)d(d) + 


dmax 
) bi(dd)S(d't) n ti(dd(d') 17 


d 
In view of the fact that the mass, and not the number, distribution is measured in practice, the popu- 
lation quantity can be changed from number to mass: e 
Let p be the density of the particles and C,, the shape factor so that the mass of a particle of 
linear dimension d is p C, d*. The massa of particles of sizes d to d+d(d) is: 
pC,d'n f(d,t) d(d) 


This same mass can be expressed in terms of the mass distribution density function: 
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(M/V) f, (d,t) d(d) 
Equation 17 can be written as: 


- | Eds) ] = — S(dyt) fa (dt) + 


*dmax 
b; (d, dy S(d)fi(d's) (d,d'y d(d”) 18 
o tú 
The quantity (d,d'): b (d,d') represents the mass of fragments in the size interval d to d+d(d) 
coming from size d' divided by the total massa of particles in size d'. In other words. it is the mass fraction of 
the breakage product from a particle of size d' ending up in size d to d+d(d). The massa breakage function 
can be defined as: 


b, (d,d') d(d)=(d'/d)” pb (dd) d(d) 19 
The final form of the population balance model becomes: 


Ce) 
e (dt) = — S(dyt) (5 (dt) + 
t 
Mimaz 
by (A dS(A e fa(d td (d?) 20 
al 1 

In order to obtain a solution to equation 20 it is necessary to specify a functional form for the selec- 
tion and breakage functions. In practice, grinding behavior is usually measured with sieve analyses, thus 
suggesting that a size-discretized form of the model would be more convenient than the size-continuous 
form. 

Define S (t) as the size-discretized selection function giving the fractional rate of breakage of parti- 
cles out of the ith size class at time t, and b,| as the size-discretized breakage function giving the fraction 
of breakage product from the jth (> ith) size interval appearing in the ith size interval. Tt is necessary 
that the size classes be small enough so that the grinding beavior is constant within the class. It has been 
found that classes interval ratios of v2 are sufficiently small for most applications. From these considera 
tions the size-discretized population balance model for batch grinding becomes: 

il 


IAW 
SD epa 4 2 bij Sj(t) AWj 21 
j 


dt 


where 4 W, = mass fraction of total material in size interval i. 
Equation 21 can be written in matrix form as: 


Saw [1-0 ]saw : 
where 
Aw: | by O So o 
aw = ja | O ás 
; é À o | 5n-1 
as AANV O Lilo Dbi cs x Bi E O. 


If S is constant and therefore independent of the size distribution in the mill which varies with time, 


then equation 22 has the analytical solution: 


Aw (1) =exp| —(1—BySt [Aw(0) 28 
” The matrix exponential defined by the infinite series: 
Es : k 
sp [—(umseJ= D[-(1-2)s: | q 
k=—o 


I-(1-B)st= 5] (1-B)S | +. 24 
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is not easily evaluated. A more convenient solution arises by making a transformation: 


AW = TAW?* 25 
where T is the matrix of eigenvectors of -(I-B) & 

() 1] 

1 a 

i—l oe ei 
Ty = s Dik 9 L; >] 

Ee 

| k=j ! 


It is necessary that S, - 5 for |-*J to evaluate T, however in the improbable case of equality, the 
appropriate values of S, and S can be made slightly different without prejudicing the simulation substan 
tially [31]. 

Substituting in equation 22 for AW: 


d E 
a (Tawe ) =—| 1—B |sTAWº or 27.a 
í E, == = de ds 


d “A 
—— AW* = T el | —B)s TAW:* 27.b 
dt — -— droe-s = —— 


The eigenvalues of —(I-B)S are distinct and equal to the S,. The right hand side of equation Z7.b 


Sa 


can accordingly be written as: 


1 
T|-(1-B)S |Tawr-nawº 28 
where À is the diagonal matrix with elements — S, 
Therefore, 
(6) 
PR e lá 
AWº(t) = - =>J(t) 29 
| O o 
e 


Back transforming to the original variable: 
AW (1)=TI(T"AW(o) 30 


When the selection functions are constant, the model is called linear. The breakage functions are not; 
and the cumulative breakage function for many systems are size normalizable, thus eliminating the neces- 
sity to separately determine the n(n-1)/2 values of the breakage function. Likewise the selection func- 
tions are often related by a power law, so that the n-1 values do not have to be separately calculated. 
Linear behavior has been verified for dry batch ball milling, but wet batch ball mill grinding and dry batch 
rod mill grinding show non-linear behavior. 

Several efficient schemes have been developed for experimentally determining the selection and brea- 
kage functions for linear batch grinding systems. In addition, the variation of these parameters with ope- 
rating conditions and dimensional factors (scale-up)) is usually of interest. These latter considerations are 
two of the weaknesses of the models as formulated, since there is no explicit indication how the selection 
and breakage functions depend on the operating conditions or mill dimensions. These relations must be 
determined empirically. 

For linear grinding systems each size class can be analyzed separately, however the accuracy of mea- 
surement decreases with particle size. This is not toomuch of a problem, since the variation of the selection 
functions with respect to particle size is linear when plotted on log-log paper. An extrapolation to the 
smaller sizes gives a good approximation to the fine-size selection functions. The parameter itself is esti- 
mated by the slope of the In Aw —vs— time linear relationship. 'This was done for dolomite in a dry labora- 
tory ball mill [4]. Figures 1 and 2 show the verification of linear behavior and power-law dependence on par- 
ticle size, respectively. 
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Figure 1 
Determination of Selection Function for 7x8 + feed 
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Figure £ 
Selection Function Dependence on alze 


The selection function for linear behavior is: 


d | 
5 = de InAW,; ( monosize particles present) 31 


E = X; X, 1 | «fa 
S= 8 E 
ud 1 o! 


Therefore, only a few of the larger size classes need to be determined, the rest being extrapolated. 
The breakage functions can be estimated as follows: consider the second (to largest) size class. From 
equation 21 we have: 
dAWa 
dt 


— Sa AWa + Da 54 AW, 33 


Remembering that initially there was only material in size class 1, it is possible to write: 


F 
% = 
do Pa Sa AWI 34 
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Therefore 


E e ada lim dAWa 
E a So 
SiAWa(lo) t>0 dt 
In general, 
SP AW, (0) [=> 0 dt 


Equation 36 implies that the cumulative breakage function B,, the fraction of material broken out 
size class 1 which ends up in size class i and smaller, can be represented by 


Ba = do lim d Fa 97 

Sp bro dt 

1 

where Fi = 1/AW(o) > AWk | the cumulative mass fraction finer than size class i. Figure 3 shows 

k==0 
for small sizes the production rate of mass is constant. This has been referred to as the zero-order production 
of fine sizes [6]. In contrast, the disappearance of the top-size is first order. The intermediate sizes at first 
increase due to arrival of particles from the largest sizes, and then decrease due to breakage of the same 
particles out of the intermediate sizes. Therefore, the disappearance and production effects must be isolated 

in order to determine the selection and breakage functions, which respectively cause these effects. 
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gero-order Production plot for 7x9 mesh feed for Dolomite 


Al the breakage functions can be determined py changing the top size to the next smaller class, for 
example 7x9 -—> 10X 14, and using equation 36. If the ratio between size classes is constant, say v2, then. 
Bij = Bij41, ren 3B 
and no 
Bnj fue bxi 39 
k—) 
is generally employed to calculate the breakage functions for the smallest size. It has been found for several 
mineral systems that the cumulative breakage function is size normalizable, thus allowing extrapolation 
to the smaller sizes. 

Several important types of batch grinding exhibit non-linear behavior. The prime examples are wet 
ball mill, rod mill, and mixture grinding. In wet ball milling the fines are suspended in the fluid, thereby 
increasing the probability of breakage of the coarser particles. In rod mills this effect is even more pro- 
nouncied since breakage occurs along the lines where the rods make contact, thus resulting in coarser par- 
ticles protecting the smaller ones. Mixture grinding is non-linear due to the effects of the components being 
ground, in which one component is softer than the other. In general, non-linear grinding behavior is cau- 
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sed by a dependence of the selection functions and possibly (but usually not) the breakage functions on the 
changing environment within the mill. In linear systems the environment also changes, but the kinetic 
parameters are independent of the changes. 

Linear systems offer several calculational advantages, among which are the existence of closed form 
analytical solutions and the availability of efficient techniques to evaluate the model parameters. To solve 
rigorously non-linear models, numerical techniques must be used. Furthermore, the implicit dependence of 
the selection function on grind time must be determined by extensive experimental measurements. Ignoring 
the non-linear nature of the system under investigation will give precise answera for the reference point 
(time) at which the model parameters were determined, but the responses decrease in accuracy as the 
region of application deviates from the reference point. 

Several implicit linearization schemes have been developed that give acceptable approximations for 
non-linear grinding systems over a small range of mill environment changes. It should be emphasized that 
the model parameters are an explicit function of the mill environment and not a function of time. The impli- 
cit relationship with time is caused by the dependence of the mill environment on time. 

A rigorous mathematical treatment of linearization of the selection function has been made by Herbst 
and Mika [6]. The breakage function is generally not considered to be dependent on the mill environment. 
The linearization is based upon a Taylor's expansion of the selection functions around the reference point 
where the selection functions were evaluated. For each selection function: 


n+ 1] f Po: s 
Si (m)Jm(t)= S (mf?) mi =: mo (Sum mim (0) mj) | bei 
In equation 40 the «*> refers to the Ag a tion; m is the vector of mass fractions for all size 
classes. 
In matrix form equation 40 becomes; 
* o * 
9 m=5mAYy (m(t) —m ) 41 
0 O 
=|) -— S(m Im 
Vi fe Um ) l ] 


Substitution of equation 41 into 21 isa 


= (m(t t))=— Sj m, +35 bij S;m; >, an (im, — my ) 42 
where E | aSstmj | 
= paia, SE e e by Vik— V 
Dik s| | Ea TER s ( po Vik— bn) 


The matrix form of equation 42 is: 


a ml pe [is |stm' em tamo 43 


where 


As formulated, equation 43 is linear and admits an analytical solution: 


m(y)= m (0) exp ( q E [esm(ay-1 Ja (1-8) S+a |m 44 


As in the linear case, the matrix exponential must be replaced by a more practical from for compu- 
tation. The similarity transformation required for the non-linear case is in general not possible to determine 
because the eigenvalues of a complete matrix like q are not calculable. The linear case had a lower trian- 
gular matrix in which the eigenvalues were calculable, in general. To facilitate this simplification it would 
be necessary to assume that the selection function in each size class depended only on the corresponding 
mass fraction, m,, and not on the entire size distribution. 'Thus 8,=S8, (m,). 


a cp =), Vj=0 1) 45 


-k . 3 Ss m, Sm 
a | PicsA qu . = e | lo - 
ne E [o (E e) ua (CB 45) 4] 


dk 
Rx — > | du 9, Si Su — BS a | 47 
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or in matrix form: 


“=[1-B]8S 8 


in which q is now a lower triangular matrix with eigenvalues & S*. Therefore, 
— i | 


-— 1 


— -1 
4 = TÊ j e dé Cen TE eee | qq * 
A fig m(0)+| T EE Us = JE e 
Normally, the values of q are determined empirically from experimental data; however, for wet ball 


milling [7] and mixture grinding [6], it has been possible to establish a relationship between q and m *. With 
the linearized solution to the batch grinding equation, a number of a non-linear systems have been simulated. 
Figure 4 shows the predictions of dry rod milling of dolomite [8]. The reference distribution was that of 1 
minute, i.e. m*=m (1). The linearized model appears to be accurate up to about 3-5 minutes: 
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Figure 4 
Dry Rod Mill Grinding of Dolomite 


Grinding applications invariably end up at the industrial level not in the batch, but in the continous 
regime. It is therefore important to extend the prediction of size reduction from batch to continuous systems. 
This is done combining mass transport through the mill with the kinetics of batch grinding. 

Population balance models of continuous milling have been developed to handle steady state or dyna- 
mic conditions in which the distribution of material along the axis of the mill is explicitly considered [9]: 
To obtain a tractable example for the purposes of this study, a simplified, vet practical, case will be briefly 
discussed. 


In many continuous grinding systems it is possible to assume that: 
1) Steady-state conditions in material transport prevail. 

2) 'The transport processa is the same for all particle sizes. 

3) The size reduction kinetics are linear. 


Under the above assumptions, the batch kineties would describe the size reduction for each particle 
size if the time of residence in the mil were known. The product size distribution can thus be considered 
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an average of the batch responses, welghted with respect to the distribution of residence times. If E (t) dt 
is the fraction of particles which spend time t to t + dt in the mill, then the continuous mill response will he: 


MOS : 
mo e] ma (E(t) de 50 


where m is the continuous system size distribution of mill product, and mg is the batch response, Combin- 
ing 50 with equation 30, we have: 


me=TJe TT m(o) 51 
(ao) 
where SG. - f exp(-S,t)E(t)dr o2.u 
Õ 
Jc, = O 52.b 


I 
Under conditions of perfect mixing in which E (t) — — exp (— t/+), + being the mean residence time, 
T 


Jeg=1/(S,2—1) 2.c 


Equation 51 and 52-a, b imply that it the amount of time that particles of a given size spend in the 
mill is accounted for and used as a weigthing factor, the continuous solution is the same as the equivalent 
batch solution summed over all possible residence times. 

This type of continuous model is called a lumped parameter model because the residence time distri- 
bution is «lumped» or averaged into one parameter, the mean residence time. In the general case, size depen- 
dent transport must be accounted for by convection and dispersion considerations. This introduces non- 
-linearities into the model, even If the particle breakage process is linear. However, in a recent computational 
study [2] it was found that the error introduced by lumped parameter approximations was on the order of 
the experimental error associated with the measurement of the corresponding product size distributions. 

In summary, although grinding is an extremely complicated process, it has been possible in the context 
of the population balance to analyze and prediet with good accuracy the behavior of both batch and conti- 
nuous syze reduction systems. In the case of mill-environment-independent kinetic parameters and material 
transport, analytical solutions are possible. These have been applied with good results to many grinding 
systems. 


FLOTATION 


The modeling of flotation on the industrial scale is a relatively recent accomplishment, the most com- 
plete models emanating from the National Institute for Metallurgy in South Africa [10]. 

In that formulation three mineral properties are considered to be the most influential on the rate of 
flotation. These are the surface activity, mineral composition, and particle size. The overall rate is also 
effected by operating properties such as the froth return rate, bubble surface area available for particle 
adherence, distribution of solids and bubble residence times. These operating properties are determined by 
the operating conditions and the mineral properties. Based on these considerations a first order rate equa- 
tion has been developed and tested in which the overall rate of flotation can be expressed as: 


E (E.G.D) = YROG(DIASL(RGD) W 58 
where r == rate at which mass enters the stream, a function of surface activity, mineral compo- 
sition, particle size. 

Y == fraction of pulp entering froth region which does not return to the pulp. 

k = constant allowing for effects other than particle size or bubble surface area. 

g (D) == function to allow for effect of particle size. 

A == total bubble surface area per unit volume of pulp. 

8 =— average fraction of total bubble surface area available for flotation. 

W == total mass of solids in the pulp. 

f, = density function for distribution of mass among properties of the solids. 


9 (D) = (2,832! *D) exp(-s/D?) (14(D/Dç)!5) 


= % DÊ +» D = max size particle floatable, 
Do = particle size where recovery is maximum, 
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For the batch cell equation 53 can be integrated to yield: 
' gas 
W (0f(k.6,Dst) = Wo) L(k.G,D;o) exp [-kgA f PSdt] 5d 


in which it is recognized that y, S are dependent on time, 
For the continuous cell a mass balance at steady state gives: 


Min Lux kG.D) — Movr fsouT(k.G.D) = VRO(DIASWI, (k,6.D) Do. 


where MM, Masp are the total mass flowrate in the feed and tails, respectively. 
If the flotation cell is assumed to be a perfect mixer, then 


L(kG,D) = four (k,6,D). 


Thus, 


| ? Mis fan (.G,D) ga 
Movur fsouT(K.,G,D) = OTDo AO db 
| 1+k76(DJASW/Mcur 
From these models all recoveries and grades can be calculated for the product streams of each cell. 
The models are complex, but computer programs available from the authors incorporate all necessary calcula- 
tions. Comparisons with industrial and pilot-plant data have been made, demonstrating the efficacy of 
the models. 


HYDROMETALLURGY 


Hydrometallurgical processes are receiving much attention recently because of the reduced energy con- 
sumption and air pollution involved in these unit operations as compared to their alternatives. The operation 
of leaching is now an important step in the processing of many metals. Almost every copper mining company 
in the world has a copper leaching operation in conjunction with their flotation process. Other metals such 
as aluminum, magnesium, lead, zinc, uranium, manganese and rhenium are products of leaching operations. 

One of the most industrially important processes is that of flotation concentrate leaching, which 
extracts in solution metal values from particles of varying properties and, in the case of continuous leaching, 
exposure to the lixiviant. The process mechanisms are multi-faceted in that they involve mass and/or heat 
transfer, disperse phase reactions, and bulk flow. 

The final section of this work will discuss the development of mathematical models for single and 
multiple particle leaching in batch and continuous systems. 


Single Particle Kinetics 


For an individual particle, the kKinetics of leaching is influenced be the geometry of the particle and 
of the mineralization. Some particles, for instance, have their mineralization confined to thin layers on solid, 
inert surfaces. In that case the surface available for reaction is unchanging. For other particles the reaction 
can occur at a sharply defined interface which moves into the particle as the reaction proceeds, 


The first type of reaction geometry in which there is a stationary reaction front can be considered to 
consist of a system or matrix of pores on the surface of or inside the particle. The lixiviant is consumed by 
mineral species C liberating metal component B. For the purpose of example it will be assumed that the par- 
ticle is spherical, 


C, (r, t) == concentration of B at distance r from center of particle at time t, 
No (r) = flux (mass/area-time) of B at position r due to diffusion. 


R, == rate of production via chemical reaction of B per unit volume, 


A mass balance for a spherical shell at radius r inside the particle gives: Rate of accumulations of B=Rate 
entering shell-Rate leaving shell+Rate of productionby chemical reaction, Thus, 


d 
47 rNg(r) — 4T(r + dr)? Ne(r + dr) ++ Re 4Tridr = a sp(rtjáTridr Vi 


From Fick's Law: 
Ns(r) = — Dg 9Co(r,t)/dr 
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where D, — effective coefficient of diffusion for component B. 


Substituting: 
Dr q dn (Et) dn (rt) a 
Da a po(r, dBi st) 58 
r or Or ot 
The flux of B arriving at the surface of the particle is equal to that going into solution: 
Up (Ra) : E 
Ny(R) = — Dy SR VM) kg | Co(k.t) — Es ot 


where R=-radius of particle, kg — mass transfer coefficient across stagnant fluid layer adjaccent to par- 


ticle surface, C, — concentration of B in the same fluid layer. 


The particular chemical reaction kinetics for the term KR, will generate a partial differential equation 
which can sometimes be solved analytically or otherwise with numerical techniques. 


For example, let R6= —k, C À 


where k — surface reaction rate constant and A is average surface area per unit volume. Tf diffusion is 
slow compared to chemical reaction and bulk transfer into the solution, then equation 58 admits the solution 


op “S = - á nú 

Ce(r.t) = Cn(r,0) =R » CP sin AM )expl — DonirtjRe) 61 
Tr n n A Ro) ; 

EE ti a 1 / CN E 

NX pl Ret) =] — at exp |! a Da nº nº t/Rº (9 
Eita 2 n* |U | 


where X (R,it) = fraction of B leached at timet. It is assumed that Cc. = O, 
If the chemical reaction is rate controlling, the solution to equation 58 is: 


Cult) =k, Galo) At  t<tmax 


s 
Catt) = C, (t max ) t > tmax 63 
Xptt) = k, Coto) At/p, t<tmax 
Xa(t) = t > tmax 


Thus the fraction leached vs time for particles of equivalent dimensions would have the forms of fig. 5. 

The case in which there is a moving reaction [ront has been intensively studied in applications to 
copper sulfide leaching. The resulting mathematical model, called the «shrinking core model», developed for 
spherical particles proposes that there is an unreacted core of minerals containing the valuable component 
which is impervious to the lixiviant. As the reaction procedes the core recedes from the particle surface 
leaving a product layer through which reactants and products travel to reach the reaction boundary or 
surrounding solution, respectively. 


To enable the development of solvable models it is generally assumed that the velocity with which 
the core shrinks is small compared to the rate of diffusion in the product layer. 'This assumption is valid for 
moderate levels of lixiviant concentration in the solution. In addition the density of the particle is consi- 
dered to be constant during leaching. 


Let a moles of lixiviant A react with c moles of mineral C to produce b moles of valuable metal B 
in ionic form. The initial radius of the particle is R; r. is the radius of the unreacted core at time t. For 
species A and B the following mass balances apply: 


Da é [, 9Ca (rx) CA (rt)= O fe 12W 
o — E meta o m— ' == = = ih F 
or dr at 64 a 
Da d à IC (rt) Cy (rt)— O ts TR 

Ei sei sa = oe 
pé dr dr At bt.b 


Equation 64 says that a steady-state exists in the product layer at each radial distance r between the core 
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o md 2 O da 


surface and the surface of the particle. This is because the core appears to be stationary relative to a 
rapidly diffusing species. 
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Figure 5 
Fraction leached vs time 


Considering that the particle mass can be calculated by: 


Mr = T r3 65 
a reaction of first order in concentration and surface area can be expressed kinetically as: 
dMr 
e k, Ca (rt) A c/a 06 
[ 


where A = 4r?, the surface area of a spherical portion of the particle of radius R or the surface area of 
a particle of radius r, k = surface reaction rate constant, c/a--ratio expressing the reactistoichiometry, 
the moles of mineral leacheed per mole of lixiviant. 
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From these considerations the core radius can be determined as a function of lixiviant concentration: 


dr. E cja ks Ca (re, L) 67 
dt To, 


In addition Fick's Law applied at the reaction and particle surface gives the following boundary conditions: 


E C: (Tu [ ) s a 
— [Já — A = ks CA (tel) BS. 
dr 
0Ca (ret b ks Ca (ret , 
a HA po 3 ks Lia (re,t) 68h 
ÔTu a 
dA (Rat) : se, e 
— DA 3R CO = ks(CA(RS) — Cs) bB.c 
0Un (Ra) e ARO 
— Do “RR is = Ke(Co(R.t) — Cn) HS. 
where C=— concentration of species in bulk solution. 
Considering component A, equation 64-a can be solved to yield: 
Cais Cat iCu pd) — Cu) 1 8) 
SA dit) as Gaste (Cas Bet) — Ca) ]) —————== de — 
nico E inte “Lrjre—r/R |— Rfre 
in which case 
2CA CiCa(rt) 7 
=D 4 Te, )=Da mate Ua 
Ore ( ) Ea IR — “| pe 
«em Ca— Calret) : 
and =p is qe SEU TO. 
O CO A ee tes 


Substituting the above expressions into equations 68-a and €68-c and combining, a solution for 
C, (rot) is obtained: 


NT - 
Cir(ret) = Cal (1 (ums) ] + (1/Da) (re-ret/R) Ata) ) l 


The recession velocity of the core can now be calculated from equation 67: 


dre erav |( E 1 ) + a 3R | ) st 
= — (Cid) Us E] ==] TO (Te = Te” ee [a 
dt ido lo E q | 
The mass fraction of C reacted at time t is: 
Xe()=1—(re/R)8 [Ê 
AN (t) 3 / rel? dr. a 
ama do E) dt 
Therefore, 
dt Roof 1 R 1/ -2/8 
—— Ma (1) = (cjajl REstiidos À emsidem IXe())—1) 1-X 9) TD 
gos [Dm fa; T/| fe +; (( ol Ls ' o(t) | 
The above equation can be a to obtain: 


1 SR Rº 


Í LXett 2,9 Ae R a 1/8 : o 
-Xe(t)| — sd [4 ils — Csstuio wa 
5) Da ke / 9DA Lo el ) [+ a [( » 4 (9) | Cat pe (a/c) TE 


JXo(0+— 


This solution is good as long as the quatity: 
à = Ca/oc(ca) Z OS 


If this condition is violated then the unsteady state solution must be used. 
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Multiple Particle Systems 


For batch systems all particles of the same size will have the same fraction leached — vs — time 
curves, Other particles will have a rather complex re lationship to the reference particles due to the size 
function exhibited by equation 76. However, the fraction leached can be transformed into a cubic equation 
in Z2=(1-X )'s and solved explicitly for X =X (R,t). Then the fraction leached of the feed size distribution 
can be calculated by considering all sizes as follows: 


"Gm xa 
Helt)= X(dofti(djd(d) 


trial 


me] 
=1 


where d=2KR, f3(d) = density function for feed size distribution. 
For continuos leaching in a well-mixed vessel, the population balance must be considered because of 


the distribution of residence times. Under the steady-state assumption the macroscopic population balance 
model reduces to: 


E Uol laT Q3] 78 
À dt p | PLS Sa 


Defining the mean residence time + as 
a mass balance on the lixiviant concentration gives: 


. q é Pa -— dimax E 7 
Cain— Ca uuT = pu Cds Yin d(d) — po Cds y d(d) 80 


ho d era Q 


Equation 78 can be solved for 1) knowing Win , and the constants associated with d(d)/dt in equation 72. 
Unless the surface reaction is rate controlling, the solution to 78 is not analytical. To arrive at this solution, 
let the feed be uniform. These assumptions are not particularly common in applications of continuous leaching 
of concentrates; however, they are not unrealistic because the uniform feed size distribution can be used as 
a basis to discretize the real system. 

d(d) ã dr. à ks CAOUT 


Thus, = 2 a e e rag ai O 81 
dt dt Bc(a/c) a 


The recession velocity of the unreacted core being constant gives as a solution assuming 

y — | ul (dmaw— Amin) Ciin al d Ed dmax 8a 
(8) 

where n/= average number of particles per unit volume in feed; dmin, dmax = minimum, maximum par- 


ticle size in feed. 


A Dix n | exp (— (dmin- d)/K Ca our + )— exp(— (dmaz — d)/K Ca our 7) | 


pras = ii 
I—K CrourT Vix/n) lexp( — dmin/KCa ouT 7) — exp(— dmaz/KCA OUT T | 

for d Ed Umin. 

ab = bin n [ 1 — exp( — (dmar — d)/RKCA OUT T | a4 
1-—KCa our + bin [mi [ exp(— dmin/KCa out 7) — exp( — dmaxz/KCas ouTT) 

for d> dmio. 

id)=4/m , td) =det(d)) f df(d) d(d) 85 

o , ale)o TI T Re 

Cain—CAOUT= pista P(KCs ouTT) 56 


b (dmax e min) 


P(X) = (duas — dBnio) X — 3 (diinaz — dêmin) Xº + 6 (dmax — Amin) ME a 
gas G(exp(— dmax/X) no exp(—dmin/X) Xá 87 


To simulate a specific case let: 


d max = 200 mm 
dmin = 50mm 
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Evo 


Po = 2.33 X 10“ moles/cc 
k, = 9.32 x 10“em/hr 


Cc, NS O & molar 


,: variable 


The values of p. and k correspond to those of chrysacolla, CuoO. SiO,. 2H,0. The lixiviant concentra- 
tion, species A being H,SO,, corresponds to a pH approximately 0.10. 
(activity coefficient of 1). For stoichiometric solids and acid n = 24,700 particles/ce. The quantity n is 
calculated from a mass balance: 


Qmax E nes 

Mr = po — d8 nrfi(dijd(d) as 
Amin b 

where M, = total mass/cc of solids in feed. Since M. is known, ie. M, = 8 X 104 moles/cc chrysacolla, and 

in the feed: 


li (dl) =— L/(dmax era min) ay 


e max 
ni = My g [ d3 fi(did(d) “O 
* “min 

The graphs of Figure 6 show how the mode] can be used to generate detailed information about 
the leaching of the mineral, In plot A the fraction leached is shown vs + and % solids (in the feed). 
This latter quantity is specified by moles of solid used per ce of solution, which in the stoichiometric case 
corresponds to & X 10-moles/cc acid. Both the number and mass size distributions in the outlet of the vessel 
are plotted in C and D as a function of +, which is equivalent to flow rate (vol/time) if the volume of the 
vessel is kept constant. The average particle size in the outlet is also show in B. It is apparent that after 
about 50 to 60 hours, the system is very close to a state of equilibrium in which the particles are remain- 
ing fixed in size; that is, not much leading is going on because the driving force of concentration gradient 
has considerably diminished. The lefthand-side of equation 86 can be solved (iteratively) to establish C,, out 
as a function of + to be plotted. (Due to lack of time this has not been done, but is not difficult to accom- 
plish). It is also interesting to note that a high % solids is not a good idea, at least for one perfect mixer, 
since the fractional conversion would have a practical limit of about 50% in the case of 20 % solids. 

With reference to experimental verification of these continuous leaching models for concentrates, it 
appears that not much detailed experimental work has been reported for analysis. In this case the theory 
has proceded the experiment. However, due to a national interest in optimum processing of copper sulfide 
concentrates in the USA, experimental programs are underway at several western universities to test the 
models. Preliminary results are encouraging; and since the models involve very few approximations that are 
unrealistic, it is expected that these considerations will lead to maximum utilization of natural resources. 
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Completado o navio Kaimei, gerador de energia a partir das ondas do mar 


Estã a progredir no Japão a pes- 
quisa sobre a produção de electri- 
cidade pelas ondas do mar. Os tes- 
tes mundiais de primeira escala 
sobre produção de energia pelas 
ondas do mar para aplicação prá- 
tica começaram em 8 de Julho, no 
Porto de Yura, Prefeitura de 
Yamagata, a bordo do navio expe- 
rimental Kaimei, que é o culminar 
de pesquisais iniciadas em 1974 pela 
Agência de Ciência e Tecnologia 
e pelo Centro de Ciência Marítima 
e Tecnologia do Japão. 

A geração de energia das ondas, 
também chamada geração de ener- 
gia absorvente das ondas, tem o 
efeito de acalmar as ondas do mar 
ao mesmo tempo que aproveita a 
sua energia para a piscicultura que 
é levada a cabo intensamente nas 
águas costeiras do Japão. 


As ondas altas são reduzidas du- 
rante a geração de electricidade 
através da absorção de energia das 
ondas. Assim, espera-se que a ge- 
ração de energia das ondas sirva 
o duplo objectivo de gerar electri- 
cidade e possibilitar que as águas 
costeiras sejam usadas efectiva- 
mente através do seu efeito cal- 
mante, 


Basicamente, o equipamento de 
geração de emergia das ondas con- 
ta com a acção das ondas para 
empurrar o ar para dentro e para 
fora das muitas câmaras de ar 
colocadas a todo o comprimento 
de um longo corpo flutuante. A 
corrente de ar atinge uma alta ve- 
locidade depois de ser forçada a 
passar atravé de bocais e põe em 
movimento turbinas de ar e gera- 
dores para produzir electricidade. 
Ao mesmo tempo, a força das on- 
das enfraquece. 

Um modelo experimental com 
T,ôm de comprimento e 1,2m de 
largura colocado num tanque de 
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água conseguiu absorver cerca de 
80% da energia das ondas como 
energia de ar, reduzindo a altura 
de grandes ondas em cerca de 2/3. 


O equipamento gerador do navio 
experimental agora completado, 
Kaimei, destina-se a gerar 200 kw 
de electricidade. As pesquisas so- 
bre o gerador começaram em 1974, 
e têm sido feitas várias inovações 
para aumentar a eficiência gera- 
dora. Já foi investido no projecto 
um total de 500 milhões de ienes 
(cerca de 2,5 milhões de dólares), 
com outros 300 milhões de ienes 
orçamentados para o corrente ano 
fiscal. 

O Kaimei, que se assemelha a 
um mnavio-tanque, tem B0m de 
comprimento e i2m de largura e 
pesa 500 toneladas. O casco con- 
tém 11 pares de câmaras de ar 
com um gerador instalado por 
cima do ponto médio de cada par. 
A água do mar move-se para cima 
e para baixo nas câmaras sem 
fundo, de acordo com a acção das 
ondas. Quando uma onda se le- 
vanta, o ar é empurrado para den- 
tro de uma das duas câmaras, e 
quando a onda cai, o ar é sugado 
para a outra câmara, A repetição 
deste processo mantém o movi- 
mento giratório da turbina de ar. 

No corrente amo fiscal, serão 
conduzidas experiências com 3 ge- 
radores instalados, respectivamen- 
te, na proa, na popa e a meia-nau, 
Sob um acordo internacional, será 
realizado no próximo ano fiscal um 
estudo comparativo, com geradores 
desenvolvidos pela Grã-Bretanha e 
Estados Unidos, a serem instala- 
dos a bordo do Kaimei. 

Eventualmente equipado com 11 
geradores, o Kaimei será capaz de 
gerar cerca de 2000 kw de electri- 
cidade. 


O uso da energia das ondas é um 


dos temas de pesquisa da Agência 
Internacional de Energia (ATE). 
Além do Japão, participam tam- 
bém neste projecto os Estados Uni- 
dos, Suécia, Canadá, Peru, Holan- 
da, Dinamarca, Bélgica e Irlanda. 

O sistema de turbina de ar do 
Japão e o método Salter Cam da 
Grã-Bretanha são considerados al- 
tamente promissores na geração de 
energia a partir das ondas. Por 
isso, são colocadas grandes espe- 
ranças na experiência a grande es- 
cala do Japão como uma grande 
contribuição no campo da utiliza- 
ção da energia das ondas. 

Uma vez postos em uso prático, 
os geradores de energia das ondas 
servirão como uma fonte de ener- 
gia para as principais ilhas, assim 
como para as mais afastadas. O 
seu efeito redutor ds ondas ajuda- 
rá também a proteger a pesca cos- 
teira e a cultura de peixes, evitar 
desastres naturais na beira-mar e 
conservar as obras de engenharia 
civil ao longo da costa. 

Particularmente digno de nota é 
o plano para fazer uso das ondas 
de inverno do Mar do Japão e usar 
a energia gerada pelas ondas para 
fazer funcionar o equipamento de 
derreter neves nas regiões de gran- 
des nevões ao norte de Honshu. 

Várias experiências a bordo do 
Kaimei continuarão até 1980, e se 
elas provarem, a geração de ener- 
gia a partir das ondas será posta 
em prática em 1981. 

Para o Japão, um pais insular 
com 13000 km de linha de costa, 
as ondas constituem a mais abun- 
dante fonte de energia ainda por 
ser utilizada. Em breve, o Japão 
será capaz de gerar 1,4 biliões de 
kw de electricidade por ano, a par- 
tir do mar que o rodeia, 


(Extraido do «Boletim Informativo 
da Embaixada do Japão», (399) Dezem- 
bro 1978). 
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Le parametrage des reserves 


RESUMO 


Distinguem-se dois aspectos da parametrização 
prática das reservas de um jazigo — a parametriza- 
cão técnica, em que se pretende obter o valor do 
teor de corte que maximiza a quantidade de metal 
para uma tonelagem fixada e uma optimização eco- 
nómica. 

Desenvolve-se a parametrização técnica tratando o 
caso da selecção livre e o caso da selecção sob cons- 
trangimentos geométricos que podem referir-se a 
jasigos estratiformes e a desmontes a céu aberto. 


Le présent texte est une synthêse de la note 
N-453 de G. MATHERON publiée au Centre de 
Morphologie Mathématique de TECOLE DES MINES 
DE PARIS em Décembre 1975. 


| — INTRODUCTION 


L'objectif théorique — et idéal — du parameétrage 
des réserves dun gisement serait de daterminer, en 
fonction des divers paramêtres qui interviennent dans 
la formule de valorisation (cours des métaux, coúts 
divers...) la meilleure exploitation techniquement 
possible, 

Formellement, WY b= (Db, «Db ) la valeur des 
paramétres économiques, il existe un ou plusieurs 
projects o(b) qui maximisent la formule de valorisa- 
tion Blw,b) à b fixé. On note y(b) le maximum. 


V(b) == B(m(b), b) = Sup=B(m,b) 


Si le projet w présente les caractéristiques techni- 
ques X, (teneur de coupure, taux de découverture, 
tonnage minerai...), on a le paramétrage: 


t(D) = X, (m(b)) avec Vb)= SupgB(m,b). 


Malheureusement, dans la pratique, il est impos- 
sible de formuler le problême de cette maniêre, pour 
2 raisons: 


— Le gisement n'est pas connu, mais seulement 
estimé. 
— Il est impensable de pouvoir optimiser à la 


D. GUIBAL * 


ABSTRACT 


Two features of the pratical ore reserve para- 
metrization of an orebody are developed — the the- 
chnical parametrization, in which the aim is to cal- 
culate the cut off grade that maxinvizes the metal 
content for a fixed tonnage and the economical opti- 
mization. 

The technical parametrizaton is emphasized in 
two cases — the free selection and the selection under 
geometric constraints in stratigraphic deposits or 
open pit operations. 


fois toutes les opérations trés complexes qui 
constituent lexploitation d'un gisement. 


L'idée essentielle qui nous guidera donc dans le 
paramétrage «pratique» des réserves sera de séparer 
le problême en phases à peu prês indépendantes, et 
de chercher, à chaque fois, des méthodes de simpli- 
fication. A titre d'exemple, déjã ancien, citons la três 
connue relation tonnage-teneur de Lasky: à partir 
d'un modéle três simple d'un porphyri-copper (teneur 
réguliêrement décroissante depuis un coeur riche), 
Lasky a cherché à étudier la variation des réserves 
récupérables en fonction de la teneur de coupure z. 


Fig. 1 


Le domaine exploité est celui limitê par Tisote- 
neur z,: 


B,. ES (x | (x) > Z,) 


Son tonnage et sa quantité de métal contenue 
sont respectivement: 


* Centre de Geostatistique de VÉcole Sup de Mines de Paris (Fontainebleau). 
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T(r)- E dx et Qu) = | Z (x) dx 


LO EA Ao dg 
O LA) 


Q (z,) 
D'oúà la teneur: m(z,) = ———. 
T iz) 

Par élimination de z, on trouve une relation 
entre Q et T et m et T. Cette derniêre peut, dans 
certaines limites, prendre la forme bien connue 
m=]—8 Log T (loi de Lasky). 

En termes d'histogrammes Fí(dz) des teneurs 
ponctuelles, Lasky utilisait la relation: 

T(dz)=T Fí(dz) ou 'T représente le tonnage ma- 
ximal possible — On en déduit les relations tonnage- 
-teneur: 


Pz)=T (1-F(z,) ) 
Q(z) = T | , ZzF (dz) 


Cet exemple — classique — de paramétrage peut 
être critiqué d'un point de vue géostatistique (Thisto- 
gramme des teneurs correspond à des carottes, non 
à des blocs; les teneurs vraies des bloes ne sont pas 
connues, mais estimées et, enfin, il existe des con- 
traintes géometriques sur la sélection qui font qu'on 
ne peut pas enlever un bloc indépendamment d'un 
autre). Tl est cependant três instructif, car il intro- 
duit lidée simple d'un paramétrage en fonction d'un 
seul paramétre technique (la teneur de coupure z ). 
Mêéme si les parameétres économiques de la formule 
de valorisation sont três nombreux, le meilleur choix 
possible du volume à exploiter en vue d'obtenir un 
tonnage T donné est celui qui maximise la quantité 
de métal Q à T fixé. 

On peut donc secinder le problême de Voptimisa- 
tion en 2 parties: 


— Un paramétrage technique: en fonction du 
paramêtre technique z , on détermine le projet 
qui maximise Q à T(z,) fixé. 

— une optimisation économique: on choisit en- 
suite z de maniêére à optimiser la formule de 
valorisation. 


C'est cette demarche — double — que Jon tentera 
de généraliser. 


| — LES FACTEURS DU PARAMETRAGE TECNIQUE : 


ls sont assez nombreux. 


a) Les caractéristiques technologiques de |'exploitation: 


Elles peuvent imposer des contraintes qui limi- 
tent le nombre de projets possibles à des contours 
autorisés. Par exemple, dans une exploitation à ciel 
ouvert, il faut respecter la pente des talus: autrement 
dit, pour pouvoir exploiter un bloc x, il faut avoir 
exploitê auparavant tous les bloes contenus dans le 
cone de sommet x et de pente celle du talus. 
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b) La sélection: 


Celle-ci peut être simple ou double. Dans ce 
dernier cas, on distinguera une sélection d'extraction 
(choix du volume V à extraire) et une sélection du 
minerai (les petits blocs v inclus dans V sont en- 
vovés, soit au stérile, soit à Jusine de traitement). 

A chaque niveau de sélection, il convient de pré- 
ciser le support (taille des bloces), et les contraintes 
yéométriques à laquelle elle est soumise. 


c) Le niveau d'information disponible : 


IH convient de distinguer soigneusement Vinfor- 
mation actuelle, disponible au moment de JVétude et 
information future, disponible au moment de la 
selection effective, qui est souvent différente de la 
premiêre et qu'il faudra anticiper en probabilite. 


dj Les critêres de sélection : 
Tis sont de la forme; 


f, (information au moment de la décision de selec- 
tion S,) >z, ou z, représente un paramétre. 

Les fonctions critiques f, doivent être telles que 
Hoptimum économique figure dans la famille de pro- 
jets définis par les f.. 

Pratiquement, on prend souvent, pour f, les esti- 


E 
mateurs Z, formés avec linformation future dispo- 
nible et les z, sont alors des teneurs de coupure. 


e) La formule de valorisation B 


Elle dépend des paraméêtres économiques et des 
caractéristiques du gisement. Souvent B sera une 
fonction de Q, (quantité de métal), T (tonnage mine- 
rai sélectionné) et V (tonnage total extrait), crois- 
sante en Q, décroissante en T et V.. 

Un cas particulicrement simple est celui oú B 
est linéaire: 


B = 20 = AV —gT 


Ill — SELECTION LIBRE 


a) Sélection libre et simple 


Dans ce cas, le gisement est découpé en bloes 
égaux v, de volume v, qu'on peut sélectionner indé- 
pendamment, il sagit de maximiser la quantité de 


métal Q =v ma A(v,), à tonnage (donc k) fixé, et 


quelle que soit la formule de valorisation B. 
Z(v,) mest pas connu: on utilise un estimateur 


Z*(v,). On impose à cet estimateur la condition de 
non biais conditionnel; 


E| 209 | 2 | == Z*(v). (condition un 


peu plus forte que la condition classique de non biais 
du krigeage, mais on peut admettre que lestimateur 
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de krigeage n'est pas trop mauvais de ce point de 
vue), 

On montre que VB, il faut utiliser comme cri- 
têre de sélection Z* (v,) = zo z paramêtre techni- 
que, est une teneur de coupure appliquée à Vestima- 
teur futur. Au stade de Iétude, Z*(v,) est inconnu; on 
devrait donc calculer pour chaque v, la quantité 
P(Z*(v;) >z | Za) (Za: information actuellement 
disponible) et estimer le tonnage total T(z) par 
Ti(z) =v 5) P(Z*(v,) > 2 | Zq). 

En fait, cette procéêdure peut être simplifiêe s il 
y a un nombre suffisant de v,: la sommation redonne 
approximativement la probabilité non conditionnelle: 


T(z) = NvP (Z* (v) > Z ) 
Pour effectuer ce calcul, il suffit de connaitre 


la loi a priori de Vestimateur futur F*(dz). 
D'oú le paramétrage cherché (en posant T =Nv) 


0 

iz) ='T, | ,  Fr(da) 
ei 

Az) =T, |, 2 F*(da) 


F*, en général inconnu, peut sobtenir à partir de 
Vhistogramme des Zy moyennant des hypothêéses de 
permanence., 

L'optimisation êconomique consiste ensuite en la 
recherche de z, valeur de z qui maximise B. 

Dans le cas três simple oú B est linéaire en 


Alv,) 
B(Z(v))=Via Zlv )-b) 


a=valeur du point métal connu 
À b=frais d'exploitation 


on obtient facilement: 


b 

Z, =— (c'est la classique rêgle marginaliste: on 
a 

exploite tout bloc qui paie ses frais en aval). 


b) Sélection libre à 2 niveaux 
On suppose à présent 2 sélections successives: 


— une sélection 8, d'extraction sur la base d'une 
information IT (pas nécessairement disponible 
au moment de Vetude) permet de choisir les 
grands panneaux V 

— une sélection S, du minerai sur la base d'une 
information I, plus riche qu'I permet de choi- 
sir dans chacun des V, précédemment sélec- 
tionné les blocs v, j à envoyver au traitement. 

|) Cas ou |, est disponible 


Au niveau de chaque panneau V, de premiêére 
sélction, on réalise le même paramétrage qu'en a); 
d'oú Vestimation (locale): 
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ão js id * 
T(z)=Vv, F, (dz | 1) 
a ds 
Pa 8] % 
Q,(2,) =V, | z F, (dz |1,) 
EA 


a E 


Ces fonctions de transfert des panneaux peuvent 
être estimées à laide du krigeage disjonctif. 

Il faut optimiser à présent la sélection S, (choix 
des V,); dans le cas simple oú la formule de valori- 
sation est linéaire: 


B=aQ-bT-eV 
V: volume extrait 


T. tonnage traité 
Q: métal contenu 


b: prix de revient du traitement 
c: coút d'extraction 
a: valeur du point métal 


On montre qu'il y a 2 paramêtres techniques (Leneurs 
de coupure): 


b 
Zz, = — valeur optimale de z à V, fixé (car cV 
est déjã dépensé) 


On effectue donc S,, pour Z et Z, selon ce critére; 
soit J E 2) Hensemble des indices j retenus, on a 
alors: 


TZat) ee > 1, (25) 
jelJ 

Q(z5321) as b> Q, (25) 
JEJ 

Vízom) = > Vi(zo) 
jJ€EJ 


Quand B m'est pas linéaire, les formules sont plus 
complexes, et on ne peut séparer le paramétrage 
technique et Voptimisation économique. 


2) Cas ou |, n'est pas disponible 


Ce cas est nettement plus difficile: il faut anti- 
ciper en probabilité les fonctions de transfert, ce 
qui conduit à TYutilisation du krigeage disjonctif. 
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IV — SELECTION SOUS CONTRAINTE 


Nous nous limiterons au cas de la sélection sim- 
ple. Un formalisme semblable, quoiqu'un peu plus 
lourd, peut être développé pour une sélection double 
(paramétrage en V et T). 

La contrainte sexprime par le fait que tout pro- 
jet doit appartenir à une famille de contours auto- 
risés. 

Ex.: ciel ouvert: un contour autorisé A doit 
respecter les pentes dequilibre (talus), donc A = 
= Ur, our, est le cône de sommet x et de pentes 


les pentes d'équilibre. 


Gisement stratiforme (toit et mur de teneur): 
dans chaque panneau P, la partie exploitge est une 
tranche horizontale de puissance > p donnée. 

Parmi tous les projets A, tels que V(A) « v, le 
meilleur est celui qui maximise la quantité meétal 
Q(A) 


Vyda, | Q(v) = Sup | Q(A) | ACBer VIA) <v ) 


Si on suppose Vinformation I qui sert au choix 
de A disponible, alors il revient au mêéme de travailler 
en teneurs vraies ou teneurs estimées (vérifiant le 
non biais conditionnel): il faut en effet maximiser 


Q'A) = | qt(x) dx 
a AM 


Malheureusement, cette maximisation n'est pas 
três facile et fait appel à une combinatoire três dé- 
veloppée. C'est pourquoi s'est développée la techni- 
que de Vanalyse convexe qui cherche à remplacer la 
courbe Q(v) par son enveloppe concave. 


G(v) 


Fig. 2 


Posons v) >0 


(04) = Sup fO(A) — GV(A)) 
A 


Soit 8, la famille des A € Bqui réalise ce sup; soient 


Ertel + - 
enfin A et A les éléments de ES. qui minimisent 
et maximisent V(A). 
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Pour VA) <v on a: QUA) — A v 4 ylA) 
et Q(v)-A v < y(A) 
avec égalité seulement si v=V(A ) pour a e Bs 
EN A 
Done yY(A) = Sup [Q(v) — À v) 
vz20 
Si on pose 
Q(v) = Inf (y(A) + av) 
1>0 


on obtient l'enveloppe concave de Q(v), comme il 
est facile de la vérifier. 


On en déduit le paramétrage des réserves en À: 


v )>0, la droite q — Av = y(h) admet deux points 
d'appui extrêmes: 


— - : + - E 
Vs et (Vo, associés à A ct A ÍJ,j,Le 
( A 1) ( EN Ec X X 


paramétrage est; 


VOO = VA) VOO) = V(A) 
q" (M)=0(4) GM) =0(A) 


dont on démontre qu'il est unique 


Donnons deux exemples: 


a) Toit et mur de teneur 


0 sagit de déterminer pour un gisement strati- 
forme les cotes a, et b, du toit et du mur en chaque 
panneau P, connaissant le krigeage des tranches 
horizontales de coté x Z*(x) 
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Fig. 4 


On optimise chacun des panneaux séparément: 
on obtient: (Z* (x) < À dês que x est au-dessus de 
a, ou en dessous de b,) 


TO) = (b-a) 


b, à 
Q,0) = ) 'Z, (ads 
da 


a 
et (4) = > T, (2) 
20) = 3, 0,0) 


Une application três détaillée de cette méthode 
à un gisement de phosphates figure dans [1). 


b) Ciel ouvert 


On a vu que les contraintes de pente sexpriment 
par le fait que dês que Von veut extraire un bloc x, 
il faut extraire tout le cone 7. On a la relation 
dordre y >x si YET. 

La famile B des contours autorisés est consti- 
tuée des ensembles stables pour 7. 

On considêre d'autre part ensemble ) des fonc- 
tions T' croissantes 
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[67 <=> (yel,=—> (9) =) ) 
On recherche les contours qui réalisent le 


sup (Q —-à Vi = v(A). 
Soit B 1a famille de ces contours, B et B”' 
A A A 
le plus petit et le plus grand élément, on associe alors 


à chaque bloc x, une fonction A définie de maniêre 
unique par: 


+ 
A(x) = Sup fà | x6B 1 
Pp (X| na) 
I est facile de voir que A€ J 


et que: B=(A>h) B =(A>A) 


On montre alors que A représente la projection 


de la teneur q(x) sur J et définit complêtement le 
paramétrage des réserves par: 


B+; = [A>3) p= [A>) 


vtj= [ NV (ds) va=] V(dx) 
Jd NDA 4 NDA 


Q*y = | AGV(d) =| n(yv(d) 


Une application importante de ce paramétrage 
est lélaboration d'un programme automatique de 
dessin de carriêre original pour lequel on pourra 
consulter [2]. 
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Encontros sobre métodos quantitativos aplicados às variáveis regionalizadas 


O Centro de Valorização de Re- 
cursos Minerais da Universidade 
Técnica de Lisboa pretende reunir, 
na 2.º quinzena de Outubro do pre- 
sente ano. um conjunto de inves- 
tigadores interessados pela aplica- 
cão de metodologias (teoria das 
funções aleatórias; processos esto- 
cásticos e estatísticas multivarian- 
tes) no tratamento de variáveis ou 
indices qualitativos detectáveis em 
sistemas naturais, 

A temática sugerida para estes 
encontros é a seguinte: 

ANALISE ESTRUTURAL — 
Análise de dados (classificação, ti- 
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pologia, taxonomia numérica, aná- 
lise factorial de correspondências, 
etc.). Espectros, covariâncias, va- 
riografia, «pattern recognition»; 

ESTIMAÇÃO E SIMULAÇÃO -— 
Estimadores lineares e não linea- 
res. Processos não estacionários. 
Simulação de variáveis correlacio- 
nadas; 

FILTRAGEM E REGULARIZA - 
ÇÃO — Atenuação de variáveis re- 
gionalizadas; 


DECISÃO E OPTIMIZAÇÃO — 


Programação dinâmica. Processos 
markovianos. Econometria; 


SISTEMAS — Aplicação da teo- 
ria dos sistemas à modelização dos 
recursos naturais. 


Os investigadores interessados 
em colaborar, deverão comunicar 
com a maior brevidade as suas cri- 
ticas e sugestões, bem como o 
tema e resumo da sua eventual 
comunicação para: 


Comissão organizadora dos en- 
contros. CVYRMUTL. IST — Pavi- 
lhão de Minas, Av. Rovisco Pais, 
1000 Lisboa, Tel. 80 12 10. 
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Dynamic programming and its aplication 
to cut-off grade problems 


RESUMO 


A determinação das sequências óptimas de teores 
de corte e ritmos de produção durante a vida de uma 
mina é um problema que tem sido particularmente 
debatido nos últimos 15 anos, 

O propósito desta comunicação é o de apresentar 
o problema no quadro da programação dinâmica: 
determinística ou estocástica, O primeiro aspecto co- 
bre os casos em que os preços e custos envolvidos 
podem ser perfeitamente fixados e o último quando 
estão sujeitos a flutuações aleatórias podendo-se an- 
tecipar o seu valor em probabilidade. 


DYNAMIC PROGRAMMING 


Dymamic programming is based on the applica- 
tion of a simple property of multistage decision 
processes. This property has been formulated by 
Bellman [1], in his principle of optimality: «An opti- 
mal policy has the property that whatever the initial 
state and initial decision are, the remaining decisions 
must constitute an optimal policy with regard to the 
state resulting from the first decision». 

At each of a given set of times a decision is to 
be chosen from a finite set of decisions. For example, 
if one of four possible decisions must be chosen at 
each decision time and the process consista of, say, 
ten such decision times, there are 4" possible different 
sequences of ten decisions each, If there is a value 
associated with each of these decisions we may be 
interested in finding the decision sequence which 
has maximum or minimum value, Such a problem 
can be solved by applying the technique of dynamic 
programming. 

The method is characterised by: 
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by PETER DOWD 


University of Leeds 


ABSTRACT 


The determination of optimal sequences of cut-off 
grades and production rates over the life of a mine 
is a problem that has received much attention in 
mining literature, especialhy during the last 15 years. 
The aim of this paper is to present the problem as a 
determinístio ar stochastic program, The former re- 
presents the case when future prices and costs can be 
perfectly predicted, and the later when they are 
subject to random changes and can be predicted in 
terms of probabilities. 


a) A function defining the total value at each 
decision time of the system, 

b) The set of permissible decisions. 

ci The present state of the system. 

d) The number of steps remaining before the 
procedure is terminated, 


Consider, for example, the simplified situation 
ilustrated in Fig. 1. A small deposit can be mined at 
either of two cut-off grades: the line of higher slope 
represents the higher cut-off grade. If the entire ope- 
ration is carried out at the lower cut-off grade, the 
deposit will be depleted after six periods, and 600 
tons of ore above the cut-off grade will be produced. 
On the other hand, the deposit will be depleted after 
two periods, if mining is carried out at the higher 
grade, and 200 tons will be produced. Between theze 
two extremes several combinations are possible, as 
shown in Fig. 1. Suppose that 100 tons of ore at the 
lower grade cam be sold at 8500, and that the same 
amount at the higher grade will fetch 81000: with a 


d21 


discount rate of 20 %, the values of each decision 
are shown in Fig. 1, 


Depletion of [ 
Ritos e De 


Quontity mined to produce 100 tons 
above cutoff grade 


| 
| 
0 


300 L00 500 500 tons 
1 À 3 É 5 6 
Period, jeors 
Higher Am 
grade grade 
Figure 1 


Decision sequence graph 


The problem is to find the route of maximal value 
between the point of departure (0) and the upper 
horizon, representing depletion of the deposit, There 
is only one route from O to A and one to B, ALA we 
can write the value 417 and at EB 833, representing 
the value of the routes to these points. We must now 
look for the routes from O to €C D and E; but, as a 
consequençe of the principle of optimality, it is only 
necessary to consider the routes from A and B that 
have already been assigned values. The route of maxi- 
mal value from O to D iz thus OBD, and at D we 
can write the value 1180. Next we look for the opti- 
mal routes from O to F, G and H and we need only 
consider those from C, D and E, and so on. 

The optimal route is thus OBDGK, i.e. to mine 
at the higher cut-off grade during the first period 
and thereafter at the lower cut-off grade. Production 
at the higher grade during the second period is traded 
off against the discounted amounts at the lower grade 
during the next three periods. It should be understood 
that this is a very elementary example, but it illustra- 
tes the basic concept of dynamic programming. 

An expression for this principle of optimisation 
can be formulated as follows: 


Let EO Sd Ho DA cre Ba DO) EE (HS) + 
á Ba (Ka, Yo) RE ram En (x, Ya! 

be the retum from producing x, tons at a grade 

Y, % in period one plus the return from producing' X, 


tons at Y. %o in period two, ete. The object is to 
maximize the retum K subject to the constraints 


n 
x,> 0 and > X, =X 
[=1 
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where x is the total ore reserve. The basic functional 
equation for this maximization is 


E: (x) =max [E (X,y)+£ (X-x)] 


Xoo NY 
n n 


where f, (x) is the benefit obtained from the optimal 
policy for n periods. Roman [15] was the first to 
apply dynamic programming to mining production 
problems but he limited his study to variable produe- 
tion rates, Noren [13] also used dynamic program- 
ming, but assumed constant average grade and cons- 
tant production rate within defined sub-regions, 


ASSUMPTIONS AND LIMITATIONS IN APPLYING 
DYNAMIC PROGRAMMING TO CUT-OFF GRADE 
PROBLEMS 


It is assumed that the distribution of the grade 
values will apply uniformly throughout the space in 
which it is defined so that, for instance, any single 
ton mined will have exactly the same grade distri- 
bution as the entire space. 

It should be noted that this is an umnrealistic 
assumption as it implies that any sample regardless 
of size or position will always have the same mean 
grade with the same variance. However, if the grades 
of a volume v are distributed lognormally within a 
portion V of a deposit then so are the volumes v' (of 
different size and shape) in the portion V' of the 
same deposit. But the lognormal parameters, median 
and variance, take different numerical values as a 
funetion of the size and shape of v and V (see «La 
theorie de DeWijs et la permanence de la lognor- 
malite» in [12]. Thus it is not the permanence of 
the distribution that is in question only the parame- 
ters of that distribution. But so long as we are only 
concerned with an «average» situation and are not 
interested in «in situ» blocks of ore the assumption 
seems valid. 

One limitation on dynamic programming solutions 
is the computer memory and time involved, A cut-off 
grade problem involving 5,000,000 decisions for exam- 
ple, takes approximately 40 minutes on an IBM 360/ 
50. However, problems can quite often be reduced 
by reducing the range of possible decisions, 

It is implicitly assumed that the entire deposit 
as defined by the grade distribution should be mined 
out, but only that proportion of material mined that 
is above the cut-off grade should be sent to the con- 
centrator. 

These assumptions are obviously more justifiable 
in the case of open pit mining than in underground 
mining. For the latter, however, their effects can 
be somewhat mitigated in that the deposit can be 
divided into a number of different distributions 
through which mining progresses, In addition, va- 
riable mining costs can be applied to any specified 
proportion of the material through which mining is 
carried out to represent the ratio of the amount of 
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ore mined to that actually sent to the surface for 
milling. 

To simplify the problem it will be assumed that 
there is no stockpilling of ore and that, at least on 
average, all ore mined within any twelvemonth period 
will also be processed during that period. It is possi- 
ble, however, to take stockpiling into account with 
the method. 


THE DEPOSIT 


As an example, a porphyry copper deposit with 
a grade distribution shown in Fig. 2 is used. This 
distribution is that of 100 ft, X 100 ft, X 100 ft, blocks 
in a 2000 ft. X 2000 ft. (horizontal plan) X 500 ft, 
(vertical) field. The mean grade of the blocks is 
1.19 %. The analysis presented here is further simpli- 
fied, in that stripping ratios and capital costs will 
not be considered, 


tons x 10º 


tons x 10º 


Relativo Cumulative 


Figure £ 


Distribution of blocks in deposit 


Concentrator operation: recovery and concentrate 
grade, 

The concentrator has been assumed to have a 
maximum throughput of 30,000 tons/day. Recovery 
has been expressed as a function of head grade and 
concentrate grade in a way felt to be typical of such 
relationships for porphyry copper deposits. In the 
example presented here, head grade may vary from 
1.19 to 1.50 %, depending on the cut-off grade, and it 
is assumed that concentrate grade will remain more 
or less stable at 25 %. 


The recovery relationship is expressed as 


Recovery, % = -—(h-1,4)? X 29.53 + 92.0 


for h< 14% 


= D2.0 for h > 14% 


where h is the head grade, and the concentrate grade 


is assumed to be fixed at 25 %. In fact, any function 
at all may be specified, including a constant recovery 


TÉCNICA 451/452 


rate. Concentrate grade may also be specified as a 
function. 


In general, the profit function cannot be ex- 
pressed as a linear function of the tonnage and crude 
ore cut-off grade but, rather, as a function of the 
tonnage and the grade of the concentrate vielded by 
the mill. Thus, the definition of recovery and con- 
centrate grade relationships is very important; these 
functions will have a significant effect not only on 
the determination of optimal schedules but also on 
the sensitivity of optimal solutions to changes in pa- 
rameters — such as costs and prices, 


ECONOMIC CONDITIONS 


Perfect competition (i.e. with the company's 
output too small to affect the market conditions) has 
been assumed in the examples presented in this paper, 
although this is not a requirement of the method. 
Provided that the information is available, prices 
can be specified as a function of production (output 
placed on the market). Prices and costs can vary 
over time, either deterministically as specified values 
or specified functions of time, or stochastically in 
terms of probability. By use of these options we can 
allow the influence of decisions at all future times 
to be felt in the current decision. Variable or constant 
taxation schedules can also be specified as part of 
the input, 


EXAMPLE: STABLE MARKET CONDITIONS 


The prices, costs and operating conditions used 
in this example are those specified im Table 1 together 
with the previously specified recovery function and 
concentrate grade. 


Table 2 shows the optimal schedule when varia- 
ble mining costs are applied to all excavations (i.e. 
the amounts shown in the column headed «Depletion 
rate of depositr): and Table 3 shows the optimal 
schedule when these costs are applied only to the 
actual amount sent to the concentrator (i.e. the 
amounts shown in the columm headed «Production 
rate»). The former case corresponds more closely to 
open-pit mining, whereas the latter is closer to under- 
ground mining — although the apportioning of varia- 
ble mining costs in underground mining would nor- 
mally lie somewhere between the two extremes, and 
would vary throughout the deposit, depending on the 
amount of material through which mining must pro- 
gress to obtain ore above the cut-off grade. 

The apportionment of the variable mining costs 
could also be interpreted in terms of selectivity. The 
application of variable mining costs to the production 
rate corresponds to perfect selectivity, whereas the 
application of these costs to the depletion rate corres- 
ponds to a more expensive selectivity in mining. 


Sds 


TABLE 1 
Data used in example 


E E es r === a er SS Et 


&875 000 per year 
$1.10 per ton 
8875 000 per year 
80.725 per ton 


Fixed mining cost 

Variable mining cost 

Fixed concentrating cost 
Variable concentrating cost 
Smelter treatment charge &185 per ton of concentrate 


Transportation charge $10 per ton of concentrate 


Gross market price for metal in concentrate $0.65 per lb 
Deductions from market price 80.05 per lb 
Deduction from concentrate grade 1.0 % 


Maximum depletion (mining) rate 36 000 tons per day 


Minimum tonnage interval for which a decision to 


mine or not can be made 1500 tons 
Maximum daily production (concentrator throughput) 30 000 tons 
Interval between cutoff grade decisions 0.025 % 
Discount rate 2o Go 
TABLE 2 


Variable mining costs applied to depletion rate of deposit 


Year Depletion rate of deposit, Production rate, Cutoff grade, % Discounted profit, 5 
tons/day tons/day 

1 33 000 30 000 0,865 17 0 384 
2 33 000 30 000 0.865 14 352 307 
3 33 000 SO 000 0.865 11 481 846 
4 31 500 SO QUO 0.815 & 155 987 
) 31 900 SO QUO 0.815 7324 790 
5 31 50H) 30 000 0.815 o 859 832 
7 31 500 30 000 0.815 4 687 865 
8 3] 50H) SO CHHO O 815 3 750 292 
y 31 500 30 000 0.815 3 UM 234 
10 30 000 SO QUO 0.640 2 364 936 
1 SO 000 SO CUM) 0.640 1 891 949 
12 50 000 30 000 0.640 1513559 
13 SO 000 30 000 0.640 1210847 
14 50 000 30 000 0.640 968 670 
Total discounted profit 585 503 506 
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TABLE 3 


Variable mining costs applied to production rate 


Year Depletion rate of deposit, 


Discounted profit, 5 


Production rate, Cutofl grade, % 
tons/day tons/day 

1 36 000 30 000 0.940 19 711 415 «+ 

2 36 000 30 000 0.940 15 769 132 

o 36 000 30 000 0.940 12 615 306 

4 56 000 30 000 0.940 10 092 245 

5 36 000 30 000 0.940 8 73 796 

6 38 000 30 000 0.940 6 459 037 

7 36 000 30 000 0.940 5 167 229 

8 36 000 30 000 0.940 4 133 783 

9 33 000 30 000 0.865 3 177 474 
10 33 000 30 000 0.865 2 941 979 
11 31 500 30 000 0.815 1 972 386 
12 30 000 30 000 0.640 1 514 042 
13-—13.8 30 000 30 000 0.640 884 373 


Total discounted profit $92 112 197 


In the schedules of Tables 2 and 3 the cut-off 
grade is lowered progressively throughout the life 
o the mine, but at different rates. The cut-off grade 
is lowered at a slower rate when variable mining 
costs are applied to the production rate than when 
they are applied to the depletion rate, The same holds 
true for the mining or depletion rate. 

Within the context of long-term stable market 
conditions it is now possible to examine the effects 
of a change in value of any of the specified variables. 
since price is generally considered to be the most 
significant variable, the effects of an increase to 
$0.75/1b and a decrease to $0.60/1b will be examined, 


When variable mining costs are applied to the 
production rate, the schedules of Table 4 are obtained, 
The mean cut-off grade is lowered to 0.859 % when 
the price is increased to $0.75/1b and raised to 0.929 % 
when the price falls to $0.60/1b. Both cut-off grade 
and depletion rate become more stable as the price 
is lowered. On the other hand, when variable mining 
costs are applied to the depletion rate a rise in price 
allows the depletion rate to increase, but, to maintain 
concentrator capacity, the cut-off grade must be 


raised as shown in the schedules of Table 5. 


TABLE 4 


Production schedules as a function of metal price: variable mining costs applied 


to production rate 


Price $0.60/1b 


Price $0.75/1b 


Depletion 
rate of Produetion Cutoff 
Year deposit, rate, grade, 
tons/day tons/day So 
1 36 000 30 000 0.940 
2 36 000 30 000 0.940 
3 36 000 30 000 0.940 
4 Su 000 30 000 0.940 
5) 36 000 30 000 0,940 
6 6 000 30 000 0.940 
T 56 000 30 000 0.940 
8 36 000 50 000 0.940 
9 S6 000 30 000 0.940 
10 56 000 30 000 0.940 
11 34 500 SO 000 0.915 
13 33 000 30 000 (1.865 
13 31.200* 30 000 0.790 


Depletion 

rate of Production Cutoff 
deposit, rate, grade, 
tons/day tons/day Em 

36 000 30 000 0.940 
36 000 50 000 0.940 
36 000 30 000 0.940 
36 000 30 000 0,940 
36 000 SU 000 0.940 
36 000 30 000 0.940 
36 000 30 000 0.940 
35 000 30 000 0.865 
d3 000 30 000 0.865 
Sl 500 30 000 0.815 
St 000 30 000 0.640 
30 000 30 000 0.640 
30 000% 30 000 0.640 


Discounted profit $47 712 448 
Mean cutoff grade 0,929 % 


Discounted profit $181 161 052 
Mean cutoff grade 0.859 % 


* Terminates at 13.3 years t Terminates at 13,9 years 
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TABLE 5 


Production schedules as a function of metal price: variable mining costs applied 


to depletion rate 


— 0000000. 
Price $0.60/1b Price 80.75/1b 


Depletion Depletion 
rate of Production Cutoff rate of Production Cutoff 
Year deposit, rate, grade, deposit, rate, grade, 
tons/day tons/day o tons/day tons/day Sa 
À 1 500 30 000 0.815 do 000 30 000 0.865 
2 31 500 30 000 0,815 do UM) 30 000 0.865 
3 31 500 St OU 0.815 da QUO 30 000 0.865 
4 31 500 SU 000 0.815 do O00 30 000 0.865 
o 31 500 SU 000 0,815 do QUO S0 000 0.865 
6 31500 30 000 0,815 do QDO Si 000 0.865 
7 30 000 30 000 0.640 SL 500 SO 000 0.815 
8 30 000 30 000 0.640 31 500 SO 000 0.815 
E 30 000 30 000 0.640 31 500 SO 000 0.815 
ló 30 000 30 000 0.640 30 000 30 000 0,640 
11 St 000 30 000 O.640 30 000 30 000 0,640 
12 SO 000 SU 000 0.640 30 000 SO 000 0.640 
13 S0 000 30 000 0.640 30 000 30 000 0.640 
14 Sl 000* 30 000 0.640 30 000% 30 000 0,640 


Discounted profit 841 376 506 
Mean cutoff grade 0.721 % 


Discounted profit $174 281 690 
Mean cutoff grade O.TT9 % 


* Terminates at 14.3 years t Terminates at 14.9 vyeara 


So far it has been assumed that all variables 95m 
wili remain stable over the life of the mine, but it is 
unrealistic to suppose that costs will not increase and 
prices will not fluctuate, The interaction of cut-off 


grades and fluctuating prices will be examined in the E 

following sections, e 85 
Ê 

VARIABLE PRICES WITH A CERTAIN FUTURE “ 
o 75) 


To examine the effect of variable prices (per- 
fectly predicted) on the operating schedule, the actual 
Canadian producers' prices between 1950 and 1970 
have been used. A constant of 41.6 cents has been 
added to each price in order to begin with a value 6s, 
of 65 cents/lb in the first year, The price variations 
are shown in Fig. 3. The results of applying these 
price variations to the previous example are shown Figure 3 
in Tables 6 and 7. 


234567 894%M0N1%4213% 151677 1819 20 
feors 


Price variations 
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Table 6 is the schedule obtained when variable 
mining costs are applied to the amounts in the rate 
of depletion column. In this case these tonnages can 
bs regarded as the actual mining rate. 

The operation begins in the first year with a 
mining rate of 31,500 tons/day and a cut-off grade 
of 0,815 %. It is then increased to 33,000 tons/day 
at 0.865 G until year 5 when, apparently in antici- 
pation of the sharp price increases of years 6 and 7, 


the production rate is lowered to 31,500 tons/day. 
Profit is actually being sacrificed in years 1 and 5 
for higher discounted profits in subsequent years. 

A more orthodox situation can be seen in Table 7, 
where variable mining costs have been applied to the 
production rate. In this case there is no immediate 
lowering of the cut-off grade with temporary price 
increases, and the sharp price increases of vears 6 
and 7 produce no change in the schedule. 


TABLE 6 


Variable prices: variable mining costs applied to depletion rate 


Depletion rate of deposit, 


Production rate, 


Year tons/dawy tons/day Cutoff grade, % Discounted profit, 3 
1 31 500 30 000 0.815 17 882 787 
2 33 000 30 000 0.865 20 791 590 
3 33 000 30 000 0.865 17 591 677 
4 33 000 30 000 0.865 15 415 108 
5 31 500 30 000 0.815 11 677 888 
6 33 000 30 000 0.865 14 140 921 . 
7 33 000 30 000 0.865 13 540 530 
8 31 500 30 000 0.815 o 890 082 
9 30 000 30 000 0.640 3 571 467 

10 30 000 30 000 0.640 3 868 323 

11 30 000 50 000 0.640 à 228 508 

12 30 000 39 000 0.640 2 408 995 

13 30 000 30 000 0.640 2154 392 

14 30 000 50 000 0.640 1633 490 

Total discounted profit $133 795 758 
TABLE 7 
Variable prices: variable mining costs applied to production rate 
Depletion rate of deposit, Production rate, 

Year tons/day tons/day Cutoff grade, % Discounted profit, 3 
1 36 000 30 000 0.940 19711 415 
2 36 000 30 000 0.940 22 289 014 
3 36 000 30 000 0.940 18 801 613 
+ 36 000 30 000 0.940 16 399 852 
5 36 000 30 000 0.940 12 545 404 
é) 36 000 30 000 0.940 14 824 670 
7 36 000 30 000 0.940 14 115 414 
8 33 000 39 000 0.865 6 142 714 
9 31 500 30 000 0.815 3 710 247 

10 31 500 39 000 0.815 4 000 869 

11 30 000 30 000 0.640 3229112 

12 30 000 30 000 0.640 2 409 478 

13 39 000 30 000 0.640 2154 778 
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Total discounted profit 8140 334 579 
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The analysis is not entirely realistic, as costs 
have been kept constant, whereas they would nor- 
mally be expected to follow rising prices, Although 
dynamic programming is capable of handling simul- 
taneous variation of all costs and prices, only price 
fluctuations have been considered — to illustrate more 
clearly the effects of such variations, Perfect predic- 
tions of prices would probably only arise from long- 
“term contracts, which would also include amounts 
to be sold at specified prices. A more realistic ana- 
lysis would be achieved if future prices were predicted 
in terms of probability, 


STOCHASTIC PROGRAMMING 


In most cases prices and costs will vary over 
time and their prediction will be subject to uncer- 
tainty, but if sufficiently accurate estimates are 
possible, they can be forecast in terms of probability. 
IL is here supposed that all costs and prices can 
assume any one of a given set of values with certain 
specified probabilities. Furthermore, it is supposed 
that at the beginning of any period each price and/or 
cost can move from one of these permissible values 
to another. It is not supposed that this change is 
deterministic; instead, it is assumed to be stochastic, 
ruled by a transition matrix P = (P,;). 

A stochastic program for prices can then be 
formulated 'as 


M 
f (x) = max > Bt Sd) Pjlm) + E, (XX) 
Des En Fa 


where p,(n) is the probability of price j in period 
n and is defined by 


M 
p,(n) = > p/(n-1) P,, j=lds M 
1=1 


and p, is the probability of a change in price from 
value i to value j. It is clear that 
0O<P;<l1 Vi) 


and 


M 
2» Pol V 


j=1 


f (x) is the optimal expected value after n periods, 
and x and y, represent optimal strategies. 


Variable prices with an uncertain future. 


Suppose, for example, that it is envisaged that 
the price will vary between $0.65 and $1.00/1b during 
the maximum expected life of the mine-say, 20 years. 
For the sake of illustration we will also assume that 
the price can only take certain discrete values bet- 
ween these two limits, as shown in Table 8, where 
the probability of each price oceurring in the first 
vear is also given. From this set of probabilities it 
can be seen that it is considered almost certain that 
a price of either $0.65/1b or 50.70/1b will prevail im 
the first year; however, there is also a slim chance 
that it will be higher. The transition probability ma- 
trix of Table 8 represents the probability that the 
price will move from one value to another in any 
given period, e.g. the first line represents the proba- 
bility that if the price in the previous period were 
$0.65/1b, in the current period it will remain at that 
value, or move to $0.70/1b or to $0.75/1b, etc. 


TABLE 8 


Transition probability matrix 


= a Ss 


Price $ 0.65 O.TO 0.75 


Probability 


at time O 0.15 0,10 
Transition U.40 0.10 
U.30 0.40 
Probability 0.10 0.30 
0,05 0.10 
0,01 0.05 
Matrix 0.01 0.035 
0.01 0.035 
0.008 0.04 


28 


0.80 0.85 0.90 0.95 1.00 
0.08 0.06 0.04 0.015 0.005 | 
0.05 0.03 0.01 0.005 0.005 
0.10 0.05 0.03 0.015 0.005 
0.40 0.10 0.03 0.01 0.01 
0.40 0.30 0.10 0.035 0.01 
0.10 0.40 0.30 0.10 0.035 
0.05 0.10 0.40 0.30 0.10 
0.05 0:10 0.40 0.30 0.10 
0.05 010 0.20 0.40 0.20 
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The resulting optimal strategy schedule when 
variable mining costs are applied to the depletion 
rate is shown in Table 9. When variable mining 
costs are applied to the production rate, the optimal 
strategy is that of Table 10, If these results are com- 
pared with those of Tables 2 and 3, it will be seen 
that the overall cut-off grade has been reduced, owing 
to the expected higher prices. By the fourth year 
there is a TT % chance that the price will have risen 


to $0.75/1b or higher, and in the fifth year this pro- 
bability has increased to 84 %. 

Only variable prices are specified here as a func- 
tion of time, but all variables may be specified as 
functional equations as, for example, in reference 
[11]. Capital and operating costs can also he consi- 
dered as functions of the cut-off grade. Higher cut- 
-0ff grades demand more selective mining and, hence, 
capital and operating costs may increase. 


TABLE 9 


Variable prices: uncertain future, variable mining costs applied to depletion rate 


Depletion rate of deposit, 


Production rate, 


Discounted expected 


Year tons/day tons/day Cutoff grade, % profit, 5 
1 da 000 SO 000 0.865 8 090 848 
2 34 000 30 000 0.865 46 472 678 
) 33 000 30 000 0.865 Si L78 142 
4 do 000 30 000 0.565 29 742 514 
o da 000 30 000 0.865 23 794 011 
Li) da 000 30 000 0.865 19 035 209 
7 31 500 SO 000 O.s15 15 062 374 
5 31 500 30 000 0.815 12 049 899 
9 31 500 30 000 0.815 9 639 919 
10) 30 000 30 000 0,640 7565 979 
11 30 000 30 000 0.640 6 052 784 
12 30 000 SU 000 0.640 4 842 227 
13 30 000 30 000 (640 o Bio Tal 
14-14.9 30 000 30 000 0.640 2 622 627 
Optimal discounted expected value of profit $276 022 992 

TABLE 10 
Variable prices: uncertain future, variable mining costs applied to production rate 
Depletion rate of deposit, Production rate, Discounted expected 

Tear tons/day tons/day Cutoff grade, % profit, & 
1 56 000 30 000 0.940 6U 267 798 
2 35 000 30 000 0.940 48 214 238 
3 so 000 30 000 0.940 8 571 391 
4 36 000 30 000 (940 30 857 112 
) 36 000 30 000 0.940 24 685 690 
(5) So 000 SD 000 0.940 19 748 552 
q Sa 000 30 000 0.865 15 464 986 
8 ds 000 30 000 0.865 12 371 989 
E SS 000 30 000 0.865 9 897 591 
10 31 500 30 000 0.815 T T64 604 
11 S1 500 30 000 0.815 6 211 683 
12 50 000 30 000 0.640 4 834 797 
13 30 000 50 000 0.640 S 867 837 
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Optimal discounted expected value of profit $282 758 267 


MINING CONSTRAINTS 


Up to this point the problem of determining 
optimum cut-off grades and production rates has been 
greatly simplified. In effect, it has been assumed 
that all ore is immediately accessible. Unfortunately, 
this is not the case. In open pit mining, for example, 
the block-grade distribution will comprise all blocks 
of ore that fall within the optimal pit contours and 
2 block can only be mined if all ore within the cone 
which has this block at its vertex has already been 
mined. To fulfill an optimal production plan, as de- 
termined by the methods outlined so far, may involve 
selecting blocks from the bottom of the ultimate pit 
during the first stages of production which would 
of course be impossible, 

Within the framework of dynamic programming 


Level 1 


Mean 114% 
Vvoronce 0063 (ho)? 


400 
300 
200 
100 
as — ai 
2 SeSo4Arcasor ss Lad ar 
un Grade,“ 
E 
[5] 
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| 
400] variance 0054 (e)! 


500 


Fig. 4 Total deposit 
divided into five levels 


350 


= — = = 


Grade, a 


it would be theoretically possible to specify all the 
constraints on selecting the blocks for every possible 
vear of production but practically this would be im- 
possible because of the enormous number of possibi- 
lities. Many simplifications, however, are possible. For 
example, during each year of production mining 
could be restricted to a maximum depth of, say, three 
blocks which would greatly reduce the number of 
constraints, Such a restriction would probably apply 
in practice anyway. This maximum vearly mining 
depth could be increased towards the bottom of the 
pit. 

As an example the distribution of figure 2 has 
been divided into the five component distribution of 
figure 4, each of which represents a depth of one 
block, i.e. a 100 ft, level. Each of these levels is to 
be mined sequentially beginning with level 1 and no 
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level can begin production until the level above it 
has been mined out. 


The resulting optimal schedule is shown in Ta- 
ble 11. Variable mining costs have been applied to 
the production rate, and prices and costs are fixed 


production, followed in the fifth year by level 3. 
Production is completed with level 5 in the thirteenth 
year. 

The programming of this problem takes up only 
slightly more storage and running time than the ori- 
ginal problem but as the number of constraints are 


TABLE 11 


Schedule for sequential mining of the five levels shown in figure 4 


Production rate, 


Year of deposit, Production rate, Cutoff Discounted Level in 
tons/day tons/day grade, % profit, 5 operation 
1 36 000 30 000 0.930 17 991 900 à É 
po 36 000 30 000 0.930 14 393 520 1 
2.4 36 000 30 000 0.930 4 605 926 1 
24-53 56 000 30 000 0.890 7 692 976 2 
4 36 000 30 000 0.890 10 257 302 2 
4.9 36 000 30 000 0.890 13895 256 2 
4,9-5 26 000 20 000 0.965 794 193 3 
6 55 C00 30 000 0.965 6 3559 588 6) 
7 36 000 30 000 0.965 o 082 866 4) 
7.3 36 000 30 000 0.965 1219888 J 
T.3-8 36 000 30 000 0.920 3 133275 4 
9 36 000 30 000 0.920 3 580 886 4 
9.7 36 000 30 000 0.920 2 005 296 4 
9,7-10 31 500 30 000 0.825 T61 048 ) 
11 31 500 SO 000 0.825 2 030 794 5) 
12 30 000 30 000 0.640 1562 397 5) 
12.6 30 000 30 000 0.640 749 950 o 


Total discounted profit $89 601 561 


at the values given in Table 1. All investment is 
assumed to have taken place before operations begin, 
and no new investment is required to open a new 
level, The probable negation of this assumption, in 
practice, will result in a slightly different schedule 
especially for the later levels, 

Production begins on level 1 at a cut-off grade 
of 0.930 % and at the maximum depletion rate of 
36,000 tons per day. These rates continue until the 
beginning of the sixth month of the third year (i.e. 
after 2.4 vears — Table 11), when level 2 comes into 
increased storage and, especially, running time will 
increase. 


Further constraints have been introduced such 
as the combination of two contiguous levels so as to 
respect open pit constraints. For example a block 
may be mined from either level 1 or level 2 but if 
mined from level 2 then the block immediately above 
it on level 1 plus the eight surrounding blocks on 
level 1 must have been previously mined (or will at 
least be simultaneously mined) as shown in figure 5. 

Each level is defined by its block-grade distribu- 
tion. Mining cannot progress to level 3 until level 1 
has been completely mined. This problem resulted 
in a 65 % increase in computing time. 


LEVEL 1 


LEVEL 2 


e 


Figure 5 
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Sol 


OTHER VARIABLES: TAXATION AND INFLATION 


The very nature of the dynamic programming 
method allows almost any variable to be taken into 
account. Various types of taxation schedules such 
as direct levies on profits and high grading penalties 
have been incorporated into the procedure, Varying 
inflation rates cam also be taken into account, 

Dynamic programming not only allows all varia- 
bles and relationships to be taken into account but 
also permits the criterion of optimisation (maximum 
present value, maximum rate of return etc.) to be 
readily changed. 


CONCLUSIONS 


Dynamic and stochastic programming are feasi- 
hle and general methods for solving cut-off grade 
and production rate problems. While the results pre- 
sented here are, to a large extent, dependent on the 
specified operating conditions of the concentrating 
and refining plants and on the definition of many 
oí the variables, it is a simple matter to alter any 
parameter either deterministically or stochastically. 

The intention in this paper was to introduce à 
general method for the solution of the cut-off grade 
problem in particular, and for mining problems in 
general, 

Much work remains to be done before a general 
theory of cut-off grades under real market conditions 
can be established; but it is felt that the methods 
presented can be used to build the foundations of 
such a theory. 
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Resumo dos Artigos publicados na «Técnica» n.º 451/452 


ama REV = — Janeiro/Povereiro/Março didi 


E 1 —— 


C. D. U. 519.24/25/27/28:622.1:550 


PEREIRA, Henrique Garcia — Introdução às variá- 
veis regionalizadas, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Fev./Mar. 1979, p. B9-96. 


Introduz-se a noção de Variável KRegionalizada no 
seu duplo aspecto aleatório e estrutural, Dão-se exem- 
plos de Variáveis Regionalizadas nos domínios geolá- 
gico e mineiro e sugere-se, em Keferências Bibliográfi- 
cas existentes no CVRMUL, outras aplicações, Inter- 
pretam-se as Variáveis Regionalizadas como realização 
de uma Função Aleatória e introduzem-se hipóteses 
restritivas. Define-se o variograma como instrumento 
básico da Geoestatística e estudam-se algumas das suas 
propriedades como revelador de estruturas. 

Indica-se o método de cálculo do variograma expe- 
rimental e listam-se alguns modelos teóricos ajustáveis. 


C. D. U. 519.25/27/28:622.1 


CORTEZ, Leopoldo Parreira — Inirodução à estima- 
ção geovestatística, «Técnica», Lisboa, LIV (457/ 
/452) Jan./Fev./Mar. 1979, p. 97-137. 


O problema da estimação geoceslatística consiste na 
inferência de estimadores de certas caracteristicas glo- 
bais ou médias de um determinado domínio (por exem- 
plo, a reserva total ou o teor médio de um jazigo 
mineral), com base na informação obtida através de 
uma amostragem parcelar, e na quantificação da pre- 
cisão dessa estimação. Esta quantificação é feita atra- 
vés da variância de estimação, e utiliza o variograma 
das variáveis regionalizadas que definem as caracteris- 
ticas a estimar, como instrumento de trabalho. 

Neste capitulo fornecem-se as funções teóricas do 
cálculo das variâncias de estimação e indicam-se alguns 
processos práticos para a realização desse cálculo. 
A Tinalizar, apresentam-se duas aplicações concretas da 
estimação geoestatistica, respectivamente à optimização 
do reconhecimento de um jazigo filoneano sub-vertical 
por galerias de nível amostradas por rouços, e ao cál- 
culo global da reserva em minério e em metal e do teor 
médio de um jazigo estratiforme amostrado por uma 
malha regular de sondagens, 


O 


C. D. U. 519.25/27/28:622 
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MUGE, Fernando Humberto — Prática de Krigagem, 
«Técnica», Lisboa, LIV (451/452) Jan./Fev./Mar. 
1979, p. 139-150. 


Introduz-se a noção de estimador linear optimizante 
— estimador de Krigagem. Constroi-se o sistema de 
Krigagem para o caso estacionário. Dão-se exemplos 
de variáveis regionalizadas em que não é válida a hipó- 
tese de estacionaridade e definem-se para estes casos 
os estimadores de Krigagem Universal. 

Dão-se exemplos de Krigagem de variáveis inter- 
correlacionadas (Cokrigagem) e constroi-se o respectivo 
sistema de Cokrigagem. 

Indica-se a estrutura básica de um programa de 
Krigagem. 


C. D. U. 519.24/25/27/28:622:550 


ROGADO, José Quintino — Uma nova parametri- 
zação de jazidas minerais, «Técnica», Lisboa, LIV 
(451/4562) Jan./Fev./Mar. 1979, p. 151-156. 


Apontam-se algumas limitações à prática habitual 
do critério de rejeicão baseado sobre pequenos blocos 
e respectivo teor tal-qual, 

É-se assim conduzido à consideração do bloco de 
beneficiamento e à construção da caracteristica optimal, 
a qual permite orientar a escolha do projecto optimal 
e, traduz, portanto a parametrização mais significativa. 

Mostra-se seguidamente que o controlo do desmonte 
é realizado após parques de homogenização e desenvol- 
ve-se a teoria básica que permite a avaliação das polfti- 
cas de constituição e retoma de tais parques. 
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Synopsis of articles published in «Técnica» n.º 451/452 


U. D. C. 519.25/27/28:622 


MUGE, Femando Humberto —- The practice of Kri- 
ging, «Técnica», Lisboa, LIV (451/452) Jan./Feb./ 
/Mar. 1979, p. 139-150. 


It is introduced the notion of best linear unbiased 
estimator — Kriging estimator. It is presented the 
Kriging system of equations for the stationary case, 

They are given examples of regionalized variables 
for the non stationary case and are defined, for these 
cases, the best linear unbiased estimator — Universal 
Kriging estimator. 

They are also given examples in which are kriged 
intercorrelationated variables (cokriging) and is pre- 
sented the respective cokriging system. 

It is shown the basic structure of a Kriging pro- 
gram. 


U. D. C. 519.24/25/27/28:622:550 


ROGADO, José Quintino —- 4 new parametrical 
description of ore bodies, «Técnica», Lisboa, LIV 
(451/452) Jan./Feb./Mar. 1979, p. 151-156. 


Some limitations are pointed out concerning the 
usual cut off criterion based on small blocks and res- 
pective crude ore grades, 

By this via the concepts of beneficiation block and 
optimal chracteristic are built up. This characteristic 
is useful for framing the design of the optimal mining 
project. Thus, it is the convenient way for the para- 
metric description of an ore body. 

Tt is also shown that the exploitation control is 
usually carried out on the output of a stockyard., The 
theoretical basis which governs the adequate stacking 
and reclaiming stockyard policies is developed. 
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U. D. C. 519.24/25/27/28:622.1:550 


PEREIRA, Henrique Garcia — Introduction to the 
regionalized variables, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Feb./Mar. 1979, p. 89-96. 


The theory of Regionalized Variables is introduced 
by its structural and random features. Examples of 
the variables in the geological and mining fields are 
presented and other applications are suggested, The 
Regionalized Variables are studied as a realization of 
Random Functions and work hypothesis are introduced. 
The variogram is definned as the basic tool of Geosta- 
tistics and some of its applications as structure depictor 
are discussed. 

The method for experimental variogram calculation 
is described and some theoretical models are adjusted, 


U. D. C. 519.25/27/28:622.1 


CORTEZ, Leopoldo Parreira — Introduction to geos- 
tatistical estimation, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Feb./Mar. 1979, p. 97-137. 


Geostatistical estimation deals with inference of es- 
timators for some global or mean characteristics of a 
certain space based on discontinuous sampling, and 
with the precision of this estimation, which is quan- 
tified by the estimation variance; the fundamental 
function in error computations is the variogram, 

Theoretical basis for calculation of estimation va- 
riances are referred, and some practical processes are 
proposed. In the last part, two examples of applications 
are given, respectively to the optimization of the sam- 
pling of a vein orebody recognized by horizontal galle- 
ries, and to the global ore reserve and mean grade 
estimations of a stratiform orebody recognized by a 
regular pattern of drill holes. 
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Resumo dos Artigos publicados na «Técnica» n.º* 451/452 


C. D. U. 519.28:622 


MARQUES, João Carvalho — Teoria e prática de 
amostragem de materiais a granel segundo o for- 
malismo de P. Gy, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Fev./Mar. 1979, p. 157-201. 


Explicam-se os mecanismos geradores dos diferen- 
tes erros de amostragem, identificam-se as variáveis de 
que dependem, estimam-se os seus dois primeiros mo- 
mentos e determinam-se as «qualidades» que uma amos- 
tragem deve possuir e como se pode, na prática, asse- 
gurar-lhas. São analisados os processos fundamentais 
de amostragem — partilha e implantação de incrementos 
— através de dois modelos matemáticos: Modelo Pro- 
babilista e Modelo de Integração, utilizando a Estuatis- 
tica Clássica e Geoestatística. 

Apresentam-se exemplos de aplicação a materiais a 
granel da indústria mineira, nomeadamente concepção 
e cálculo de instalações automáticas de amostragem. 

Tem, contudo, maior generalidade, sendo aplicável 
a materiais a granel de outras indústrias, nomeada- 
mente metalúrgica, química, farmacêutica, agricola e 
alimentar, do cimento, do vidro, da cerâmica e do pa- 
pel, a resíduos urbanos e industriais e ainda a mate- 
riais contínuos, sólidos ou líquidos, como chapas, 
barras, lingotes e produtos metalúrgicos fundidos em 
deslocamento, 


C. D. U, 519.24/25/28:622.1:550 


ROGADO, José Quintino — Representatividade da 
amostragem de jazidas minerais, «Técnica», Lisboa, 
LIV (451/452) Jan./Fev./Mar. 1979, p. 203-208. 


São analisados aspeclos da representatividade da 
amostragem de jazidas minerais destinada a ensaios 
mineralúrgicos e metalúrgicos, 

é definida população de referência e apresentada 
metodologia para a reproduzir por amostragem e apre- 
senta-se um exemplo real extraído do estudo de exten- 
sa bacia carbonifera Brasileira. 


Ano LIV — Janeiro Fevereiro Março 1979 


C. D. U. 519.28:622.1/2/7:550.4/7:553 


PEREIRA, Henrique Garcia — Introdução à simula- 
ção, «Técnica», Lisboa, LIV ([451/452) Jan./Fev./ 
!Mar. 1979, p. 209-224. 


Os princípios da simulação estocástica são discuti- 
dos na sua generalidade e o Método de Monte Carlo é 
descrito brevemente, 

Aplicações da Simulação Estocástica nos diferentes 
domínios das Ciências da Terra são referidas no texto 
e bibliografia existente no CVRMUL é indicada para 
cada caso específico. 

Os domínios considerados são a Geologia e Prus- 
pecção, Exploração e Tratamento de Minérios e fina- 
liza-se com a Análise econômica global de um projecto 
mineiro, 


C. D. U. 519.28:622 


DOWD, Peter — Simulação condicional de variá- 
veis regionalizadas, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Fev./Mar. 1979, p. 225-234. 


A teoria da simulação condicional de variáveis re- 
gionalizadas é introduzida e as propriedades da variável 
simulada são descritas. 

O método das bandas rotativas para a simulação 
não condicional a três dimensões é apresentado, bem 
como a simulação a uma dimensão obedecendo a uma 
covariância dada e a um histograma experimental. 

Os problemas de anisotropia, efeito de pepita e 
deriva são introduzidos na teoria e estudadas as mo- 
dificações a realizar para que esses aspectos estruturais 
sejam tomados em conta. 

Alguns exemplos são apresentados para ilustrar as 
aplicações da teoria, 
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PEREIRA, Henrique Garcia — Introduction to simu- 
lation, «Técnica», Lisboa, LIV (451/452) Jan./Feb./ 
!Mar. 19798, p. 2095-224. 


The principles of stochastic simulation are briefly 
discussed and Monte Carlo method is described. 

Applications of stochastic simulation in some fields 
of the Earth Sciences are referred and references are 
given for each specific case, 

The main fields treated are geology and explora- 
tion, mine planning and ore dressing. The global eco- 
nomic gtudy of a minmg venture is also introduced in 
a case story. 


U. D. C. 519.28:622 


DOWD, Peter — The conditional simulation of re- 
gionalized variables, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Feb./Mar. 1979, p. 225-234. 


The conditional simulation theory of Regionalized 
Variables is introduced and the basic properties of the 
simulated variable are described. 

The turning bands method for non-conditional si- 
mulation in a three-dimensional space is presented as 
well as the one-dimensional simulation generated by a 
given covariance function and experimental histogram. 

Problems of anisotropy, nugget effect and drift are 
introduced in the theory and modifications for taking 
these structural features into account are done accor- 
dingly. 

Some case stories are presented for application of 
the theory. 
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LIV — January February March 1979 


U. D. C. 519.28:622 


MARQUES, João Carvalho — Particulate material 
sampling theory and practice by P. Gy formalism, 
«Técnica», Lisboa, LIV (451/452) Jan./Feb./Mar. 
1979, p. 157-201. 


The paper explains all types of sampling errors 
and estimates their first and second moments. It defi- 
nes the «qualities» and the rules to respect in sampling 
procedure. 

The main sampling procedures are studied through 
two mathematical models: Probabilistic and Integration 
REGE, developed by Classical Statistics and Geosta- 
tistics. 

Applications to bulk materials of mining industry 
namely predesign of automatic sampling plants. 

The models can be applied to bulk materials of 
many other industries such as Metalurgical, Chemical, 
Food, Agriculture, Concrete, Glass, Ceramices and Pa- 
per, to industrial and city wastes and even to continuos 
solid and liquid materials such as ingots, rods, bars, 
steel plates and fused metalurgical products in conti- 
nuos mooving. 


U. D. C. 519.24/25/28:622.1:550 


ROGADO, José Quintino — Representativity of ore 
bodies sampling, «Técnica», Lisboa, LIV (451/452) 
Jan./Feb./Mar. 1979, p. 203-208. 


Comments are developed on the representativity of 
ore bodies sampling for obtaining data conceming the 
behaviour on concentration and metallurgical processes. 

The reference population iz defined and some me- 
thods for reproducing it are presented, | 

An example taken from a study carried out in a 
Brazilian coal basin is summarized. 
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C. D. U. 338.5:519.28:622 


SANTOS, Heitor Simões — Avaliação económica de 
projectos mineiros por simulação, «Técnica», Lis- 
boa, LIV (451/452) Jan./Fev./Mar. 1979, p. 235- 
-260. 


Dado o carácter prático deste trabalho, e estando 
ele integrado no Seminário de Geomatemática, onde 
outros trabalhos apresentam os conceitos teóricos da 
simulação, estes aspectos não são incluidos aqui. 

Por outro lado, sendo este processo de simulação 
uma aplicação à Análise Econômica de Projectos admi- 
te-se que: 


— São conhecidos os conceitos necessários ao 
Discounted-Cash-Flow 

— Um projecto/exemplo está disponivel e já foi 
explicado. 

— São conhecidos os princípios básicos da Tn- 
formática. 


Para satisfazer os dois primeiros pontos o leitor 
deve remeter-se ao anexo, onde de uma maneira su- 
cinta se relembram os conceitos, se introduz o projecto/ 
fexemplo, e se analisa o efeito de variações discretas 
(Análise de Sensibilidade). 

Neste trabalho propriamente dito, expõe-se a estru- 
tura do processo de simulação, apresentando o projecto/ 
(exemplo na óptica dos conceitos probabilísticos (Aná- 
lise do Risco). 


C. D. U. 519.28:622.271 


MUGE, Fernando Humberto — Determinação do 
contorno final dum desmonte a céu aberto (ullti- 
mate Pit»), «Técnica», Lisboa, LIV (451-452) Jan./ 
/Fev./Mar. 1979, p. 261-270. 


Formula-se na gua generalidade o problema da de- 
terminação do contorno fTinal dum desmonte a céu 
aberto. 

Descreve-se o algoritmo de Lerchs e Grossman com 
base em noções da teoria dos grafos, 

Introduz-se o método dos cones deslizantes aplicado 
ao algoritmo de Lerchs e Grossmann e dá-se um exem- 
plo de aplicação, Descreve-se a estrutura básica de um 
programa que utiliza este método. 


C. D. U. 622.2:519.28 


COSTA, J. Beja e VALENTE, Jorge — Programação 
de desmontes em lavra subterrânea, «Técnica», Lis- 
boa, LIV (451/452) Jan./Fev./Mar. 1979, 271- 
-291. 


À exploração de minério de ferro da «Companhis 
Mineira do Lobito, 8. A. R. LL.» (Cassinga-Angola; 
ox<10º ton/ano de concentrados) foi programada, duran- 
te 6 anos, com sucesso na estacionarização de parâme- 
tros, tais como o número de máquinas em operação nos 
desmontes, o teor dos concentrados e o respectivo ren- 
dimento em peso, assim como na repartição de finos 
(1/5 mm) nesses concentrados. Nesta programação de 
desmonte introduziram-se ainda alguns melhoramentos, 
tais como a optimização da recuperação em peso, com a 
introdução de um algoritmo computorizado e mantendo 
simultaneamente o teor em ferro dos concentrados num 
certo valor pré-determinado; um outro melhoramento 
foi o estabelecimento de uma sequência de desmonte 
que define as combinações de blocos a serem desmon- 
tados simultaneamente ou sequencialmente, de forma a 
se assegurar as condições descritas anteriormente, 

Este método de programação é agora adaptado ao 
caso particular das minas subterrâneas e é demons- 
trada a forma de relacionar o cálculo dos diferentes 
tipos de equipamento destas minas pela primeira vez, 
segundo supomos. 

No final descreve-se um exemplo prático de aplica- 
cão do método: este exemplo refere-se à Mina do 
Moinho (Aljustrel) de «Firites Alentejanas, S.A R.L.». 


C. D. U. 622.73/765 


OBLAD, Alex — Métodos estocásticos e determinis- 
ticos de análise e simulação em processos da meta- 
lurgia extractiva, «Técnica», Lisboa, LIV (451/4552) 
Jan./Fev./Mar. 1979, p. 293-313. 


Neste artigo apresentam-se métodos correntemente 
usados na construção e análise de modelos empíricos 
pela via estocástica clássica e também recorrendo a 
modelos deterministicos recentemente criados. Nesta úl- 
tima óptica e usando o modelo de balanço generalizado 
em sistemas de partículas, tratam-se aplicações à moa- 
gem discontinua/continua, à flutuação e à lixiviação. 
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COSTA, J. Beja e VALENTE, Jorge — The trans- 
portation system planning in underground mining, 
«Técnica», Lisboa, LIV (451/452) dJan./Feb./Mar. 
1979, 271-291. 


The exploitation of the iron ores of «Companhia 
Mineira do Lobito, S. A. R. LL.» (Cussinga-Angola; 
xi PPY of concentrates) is based upon a sophisti- 
cated planning schedule which, for six vears, has been 
successful in Keeping several pertinent parameter 
steady, such as the number of operating machines in 
the mine, the iron content and the weight recovery ol 
lhe concentrates, as well as the proportion of fines 
l/) mm) in such concentrates. In this planning sche- 
dule some improvementa was introduced, such as the 
optimization of the obtained weight recovery, with the 
introdution of a computerized optimization algorithm, 
and simultianeously keeping the iron content of lhe 
concentrates at a predetermined value; another impro- 
vement was the establishment of working criteria for 
the mine which define the combination of blocks to be 
simultaneouslyv or sequentially treated in order to en- 
sure the conditions described, 

AW this method is now, in the present text, adapta- 
ted to the particular case of the underground mines 
and the way to relationate the calculation of the num- 
ber of the diferent types of equipment used in that 
mines is shown, we suppose, for the first time. 

In the end, a pratical example of the method is 
also described. this example concerns to the Mina do 
Moinho (Aljustrel) of the Company «Pirites Alenteja- 
nas, 3. A. R. Lx. 


U. D. C. 622.73/765 


OBLAD, Alex — Process analysis and simulation in 
extractive metallurgy by stochastic and determinis- 
tic methods, «Técnica», Lisboa, LIV (451/452) 
Jan./Feb./Mar. 1979, p. 293-313. 


This paper deals wilh current methods used in 
construction and analyeis of empirical models by the 
classical stochastic approach and with recent deve- 
Inpments in deterministic models. These latter include 
the applications of a generalized population balance 
model for particulate systems to the operations of 
batech/continuous grinding, flotation, and concentrate 
teaching. 


TECNICA 451/452 


U. D. C. 338.5:519.28:622 


SANTOS, Heitor Simões — Mining projects under 
risk analysis concepts, «Técnica», Lisboa, LIV (451/ 
/452) Jan./Feb./Mar. 1979, p. 235-260. 


As this paper intends to be a practical one as well 
as a part of the present semminar, where theoretical 
concepts on simulation were already introduced, the 
former part was skipped. 

On the other hand, being «This Simulation Pro- 
cess» concerned with the Economic Analysis of the 
Projects three points are assumed: 


— Modern concepts of Econ, An, of Proj, (DCF) 
are known. 

— À project/example is available and alrealy 
explained. 

— Basic concepts on the computer science are 
required. 


Excluding the last point, the reader is reverted to 
the Annex, where a short introduction on DCF is 
available (as reminder) as well as the Project/Example 
being introduced in order to clear the structure and 
the concepts explored on discrete variations (Sensiti- 
tivity Analysis). 

Therefore, the present paper iz concerned with the 
structure of a simulation process, presenting The Pro- 
ject/Example under the probabilistic concepts (Risk 
Analysis). 


U. D. C. 519.28:622.271 


MUGE, Fernando Humberto — Determination of an 
optima! Ultimate Pit contour, «Técnica», Lisboa, LIV 
(451/4552) Jan./Feb./Mar. 1979, p. 261-270. 


The problem of the determination of an ultimate 
pit contour in order to maximize some objective fonc- 
tion is formulated in all ita generality. 

It is described the Lerchs and Grossmann alghoritm 
based in concepts of the graphs theory. 

IH is introduced the moving cone method applied 
to the Lerchs and Grossmann alghoritm and it's given 
an example of application. It is described the basic 
structure of a computer program which performs the 
above mentioned method, 
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GUIBAL, D.-—-A4 parametrização das reservas, 
«Técnica», Lisboa, LIV (451/452) Jan./Feb./Mar. 
1979, p. 315-319, 


Distinguem-se dois aspectos da parametrização prá- 
tica das reservas de um jazigo — a parametrização Léc- 
nica, em que se pretende obter o valor do teor de corte 
que maximiza a quantidade de metal para uma tonela- 
gem fixada e uma optimização económica. 

Desenvolve-se a parameirização técnica tratando o 
caso da selecção livre e o caso da selecção sob cons- 
trangimentos geométricos que podem referir-se a jazi- 
gos estraliformes e a desmôntes a céu aberto, 


C. D. U. 519.28:338.5:622 


DOWD, Peter — Aplicações da programação diná- 
mica à determinação de teores de corte, «Técnica», 
Lisboa, LIV (451/452) Jan./Feb./Mar. 1979, p. 
321-332. 


A determinação das sequências óptimas de teores 
de corte e ritmos de produção durante a vida de uma 
mina é um problema que tem sido particularmente 
debatido nos últimos 15 anos. 

O propósito desta comunicação é o de apresentar 
o problema no quadro da programação dinâmica: de- 
terministica ou estocástica. O primeiro aspecto cobre os 
casos em que os preços e custos envolvidos podem ser 
perfeitamente fixados e o último quando estão sujeitos 
à flutuações aleatórias podendo-se antecipar o seu va- 
lor em probabilidade. 
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| DOWD, Peter — Dynamic programming and its 
application to cut-off grade problems, «Técnicas, 
Lisboa, LIV (451-452) Jan./Fev./Mar. 1979, p. 
321-332. 


The determination of optimal sequences of cut-off 
vrades and production rates over the life of a mine is 
a problem that has received much attention in mining 
literature, especially during the last 15 years. The aim 
ow this paper is to present the problem as a determi- 
nístic ar stochastic program. The former represents the 
case when future prices and costs can be perfecily 
predicted, and the later when they are subject to ran- 
dom changes and can be predicted in terms of proba- 
hilities. 


TECNICA 451/452 


U. D. C. 519.24/25/27/28:622:550 


GUIBAL, D. — Le parametrage des reserves, «Téc- 
nica», Lisboa, LIV (451-452) Jan./Fev./Mar. 1979, 
p. 315-319, 


Two features of the pratical ore reserve parame- 
trization of an orebody are developed — the thechnical 
parametrization, in which the aim is to calculate the 
cut off grade that maximizes the metal content for a 
fixed tonnage and the economica) optimization. 

The technical parametrization is emphasized in two 
cases — the free selection and the selection under geo- 
metric constraints in stratigraphic deposits or open pit 
operations. 
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